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Abstract 

We compute the non-MHV one-loop seven-gluon amplitudes in M = 4 super- Yang-Mills the- 
ory, which contain three negative-helicity gluons and four positive- helicity gluons. There are 

four independent color-ordered amplitudes, ( h+++), ( 1 ( H h+) and 

( — | 1 h+). The MHV amplitudes containing two negative-helicity and five positive-helicity 

gluons were computed previously, so all independent one-loop seven-gluon helicity amplitudes are 
now known for this theory. We present partial information about an infinite sequence of next-to- 
MHV one-loop helicity amplitudes, with three negative-helicity and n — 3 positive-helicity gluons, 
and the color ordering ( h+ • • • ++); we give a new coefficient of one class of integral func- 
tions entering this amplitude. We discuss the twistor-space properties of the box-integral-function 
coefficients in the amplitudes, which are quite simple and suggestive. 

PACS numbers: ll.15.Bt, ll.25.Db, 11.25.Tq, 11.55.Bq, 12.38.Bx 



2 



I. INTRODUCTION 



Gauge theories play a central role in modern theoretical physics. They form the backbone 
of the Standard Model of particle interactions. They also play an increasingly important 
role in our understanding of string theories via the AdS/CFT correspondence. That corre- 
spondence is a strong-weak coupling duality between type IIB string theory on an AdS 5 x S 5 
background and four-dimensional Af = 4 supersymmetric gauge theory. In a recent paper, 
Witten [1] proposed a weak-weak coupling duality between Af = 4 supersymmetric gauge 
theory and the topological open-string B model on CP 3 ' 4 using D-instanton contributions. 
This proposal generalizes Nair's earlier description [2] of the simplest gauge-theory ampli- 
tudes as correlation functions on CP 1 . Berkovits [3, 4], Neitzke and Vafa [5], and Siegel [6] 
have given alternative descriptions of such a possible topological dual to the Af = 4 the- 
ory. Both of these dualities motivate the computation of perturbative amplitudes in gauge 
theory. 

Indeed, the availability of an extensive literature of explicit results for both tree-level 
and loop amplitudes in gauge theories has provided an important stimulus to and guide in 
building the topological constructions. Most notable are the series of helicity amplitudes 
with arbitrarily many external gluons, the maximally helicity violating (MHV) amplitudes 
at tree level [7, 8], and at one loop in supersymmetric theories [9, 10]. These amplitudes 
display a remarkable simplicity. The entire amplitude at a given multiplicity is much simpler 
than the typical Feynman diagram, chosen from the great number that are assembled into 
the amplitude in a traditional approach. There are indications that this simplicity continues 
to higher loops [11, 12, 13]. 

The forthcoming generation of experiments probing beyond the Standard Model at the 
LHC provides another important motivation for computing amplitudes in perturbative gauge 
theory. Precision measurements at SLC and LEP have proven to be a powerful means of 
advancing our understanding of the Standard Model. A drive towards precision physics at 
hadron colliders will require a corresponding theoretical effort in precision calculations in 
QCD, to higher loops and multiplicities. While a traditional field-theoretic point of view 
would view the pure glue contributions to QCD amplitudes as a small part of the Af = 4 
result, it is more natural to view matters the other way around. The Af = 4 results are much 
simpler than the QCD ones and can be regarded as a building block for the latter [9]. For 
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example, the gluon (or quark) circulating in the loop in a one-loop n-gluon QCD amplitude 
can be decomposed into a linear combination of an Af = 4 super-multiplet (absent in the 
quark case), an Af = 1 chiral super-multiplet, and a scalar in the loop. Each piece has 
a different analytic structure, so it makes sense to separate the components. The Af = 4 
amplitudes can always be written as a sum over a more limited class of integral functions than 
the other components, as they require only scalar box integrals. These four-point integrals 
have no loop momenta in the numerator of the integral [9]. The remaining computational 
problem is to determine the coefficient with which each scalar box integral appears in the 
amplitude. This simplicity can be traced back to the much-improved ultraviolet behavior of 
the Af = 4 theory. At the same time, the Af = 4 component captures the leading infrared 
singularities present in the QCD amplitude. 

Calculations in QCD and in Af = 4 super- Yang-Mills theory are closely related to each 
other. In general, an Af = 4 amplitude may be extracted from a corresponding QCD 
amplitude, if the number of gluonic and fermionic states is tracked in the computation. The 
conversion of a QCD amplitude to an Af = 4 super- Yang-Mills amplitude is accomplished 
by assigning 8 spin degrees of freedom to gluons circulating in the loops, and modifying 
the multiplicity of the fermions to also carry a total of 8 spin degrees of freedom. The 
color algebra also needs to be modified to account for the gluinos being in the adjoint 
color representation, instead of the quarks' fundamental color representation. With these 
modifications, the theory has the number of degrees of freedom of ten-dimensional Af = 1 
super- Yang-Mills theory compactified to four dimensions, which is another name for Af = 4 
super- Yang-Mills theory. This simple conversion holds at any loop order. As a non-trivial 
example, after carrying out this conversion on the two-loop four-gluon amplitudes given in 
ref. [14], they agree perfectly [12] with the previously determined Af = 4 amplitudes [11]. The 
analytical results presented here can therefore serve as benchmarks for testing algorithms [15] 
for direct numerical evaluation of high-multiplicity one-loop amplitudes in QCD. 

The experimentally-driven computations are technically complicated, so they require 
more powerful methods than textbook ones. In the past decades, a number of new 
approaches have been developed to cope with this complexity, including helicity meth- 
ods [16], color decompositions [17, 18, 19, 20], recursion relations [8], supersymme- 
try Ward identities [21], ideas based on string theory [22, 23], and the unitarity-based 
method [9, 10, 24, 25, 26]. The latter technique has been applied to numerous calculations, 
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most recently the two-loop calculation of all helicity amplitudes for gluon-gluon scatter- 
ing [27] and the two-loop splitting amplitudes g — > gg [26]. (The latter agree with the 
explicit collinear limits of certain two-loop amplitudes [26, 28].) 

Forefront calculations still tax available technologies to their utmost, which motivates 
the continuing development of new techniques. The twistor-space representation, intro- 
duced along with the topological string theories mentioned above, shows great promise as a 
source of new methods for performing cutting-edge calculations. Indeed, Cachazo, Svrcek, 
and Witten [23] have used observations about the twistor-space structure of amplitudes to 
formulate a new and simple set of rules based on MHV vertices (off-shell continuations of the 
Parke- Taylor MHV amplitudes) which can be used to compute all tree-level amplitudes in 
unbroken gauge theory. Their rules extend to processes with external fermions and scalars 
as well as gluons [29], and lead to explicit results for next-to-MHV amplitudes [23, 29, 30] 
and checks of 'googly' amplitudes (with two positive-helicity, and the rest negative-helicity, 
gluons) [31, 32, 33]. The rules can be implemented either analytically or numerically. The 
computational efficiency of these rules may be further enhanced with a recursive rearrange- 
ment [34]. 

A critical question is how to extend these ideas to loop calculations. The challenge is to 
reduce loop calculations to a purely algebraic problem of polynomial complexity, avoiding 
the exponential explosion in the complexity of intermediate stages that would be encoun- 
tered with brute force tensor reduction and integration methods. A direct topological-string 
approach appears to lead to complications, mixing in non-unitary states from conformal 
supergravity [35]. However, in an important step, Brandhuber, Spence, and Travaglini [36] 
have shown that one can compute one-loop amplitudes using exactly the same MHV ver- 
tices that work at tree level. They showed explicitly how to compute the infinite sequence of 
one-loop MHV amplitudes in Af = 4 super- Yang-Mills theory. Many technical steps in their 
computation parallel the original cut-based method [9]; however, the new computation is 
quite different conceptually, and probes the off-shell continuation of the MHV vertices. The 
BST calculation also reaffirms the basic simplicity of the twistor-space structure of one-loop 
amplitudes. Taking into account [37] a 'holomorphic anomaly' brings an earlier investigation 
of this structure [38] into agreement with the picture emerging from the BST calculation. 
A delta function in the holomorphic anomaly for a unitarity cut completely freezes the in- 
tegration over the intermediate phase space variables [39, 40], making it simple to evaluate 
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the action of the differential operators for twistor space co-linearity or co-planarity on the 
cuts of one-loop amplitudes. The anomaly has already been applied by Cachazo [40], in 
conjunction with unitarity, to derive algebraic equations relating coefficients of scalar box 
integrals to the rational functions making up a cut integrand. 

In this paper we present new results for amplitudes in the J\f = 4 gauge theory. We will 
give compact formulae for all four helicity configurations required for the seven-point next- 
to-MHV (NMHV) amplitude. These amplitudes were computed by evaluating unitarity 
cuts in various channels, the same method used previously to obtain the all-multiplicity 
MHV helicity amplitudes [9] and the six-point NMHV helicity amplitudes [10]. In this case, 
however, the unitarity cuts are more complicated. We have reduced them to a standard set 
of (cut) scalar box integrals using integral reduction methods implemented on the computer. 
The coefficients of the scalar boxes in the full amplitude are equal to the coefficients of the 
cut scalar boxes that have a cut in the channel considered. The resulting expressions for 
the coefficients are analytical, but quite large. However, remarkably simple expressions 
exist for these quantities. The simple expressions can be found by considering the analytic 
behavior of the coefficients in various limiting regions of the seven-point phase space, and 
using that information to build ansatze for the coefficients. The ansatze can then be checked 
numerically with high precision against the large expressions at random kinematic points. 

There are a large number of box coefficients to determine (92, after invoking some re- 
flection symmetries). However, their structure can be fit into general patterns, which are 
simple to describe in twistor-space language. The twistor-space structure of the full ampli- 
tude, including the box integrals multiplying these coefficients, is less transparent because 
of the holomorphic anomaly which affects the integrals [37] . On the other hand, the twistor- 
space properties of the coefficients square well with the recent application of the holomorphic 
anomaly to computation of Af = 4 box coefficients [40] . 

We also give an all-n formula for the coefficient of a particular class of (three- 
mass) box integrals in the adjacent-minus NMHV amplitude, the helicity configuration 

( h+ • — h+). We study the twistor-space structure of this part of the amplitude and 

find that it is consistent with expectations of simplicity emerging from the BST calculation 
and the holomorphic anomaly. These results should provide a useful guide to the analytic 
structure likely to emerge in other scalar-box coefficients in this amplitude and in other 
amplitudes. 
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As this paper was being completed, ref. [41] appeared, in which one of the four seven- 
point helicity amplitudes presented here, ( h+++), is computed via the holomorphic 

anomaly and unitarity. Although some of the expressions obtained for the box coefficients 
are more complicated than ours, we have compared them all numerically and find complete 
agreement. 

This paper is organized as follows. In section II we review the color and helicity de- 
compositions of tree-level and one-loop gauge amplitudes, as well as the structure of the 
MHV amplitudes at these orders. In section III we describe the application of the unitarity- 
based technique to the seven-point NMHV computation, including an outline of the integral 
reduction approach and its implications for the denominators of different classes of box 
coefficients. We also sketch how the coefficients were simplified. Section IV details the 
results for the seven-point amplitudes by listing the independent box integral coefficients. 
Section V describes consistency checks that were applied to the results. These checks in- 
clude the behavior as two gluons % and j become collinear, i.e. as the kinematic invariant 
$ij = (ki + kj) 2 = 2ki ■ kj — > 0. They also include the behavior as multi-particle invariants 
vanish; that is, = (fcj + kj + k k ) 2 — > 0. In section VI we give one of the box integral 

coefficients in the all-n NMHV helicity amplitude ( H h+). Section VII analyzes 

the twistor-space properties of the seven-point box coefficients, and of the term in the all-n 
NMHV helicity amplitude from section VI. In section VIII we present our conclusions. 
There are two appendices. Appendix A collects the dimensionally-regulated scalar box in- 
tegral functions appearing in the M = 4 amplitudes. Appendix B describes an example of 
how one of the seven-point helicity amplitudes can be factorized onto a multi-particle pole. 

II. REVIEW OF PREVIOUS AMPLITUDE RESULTS 

We now briefly summarize results for previously computed amplitudes in H — 4 super- 
Yang-Mills theory. At tree level it is convenient to use the color decomposition [17, 18] of 
amplitudes 

AT({h,Ka l })= J2 TV(rM 1 )...r^))A t I1 >(l Al ) ... ) ^)) ) (1) 

(T£Sn/Z n 

where S n /Z n is the group of non-cyclic permutations on n symbols, and j Xj denotes the j-th 
momentum and helicity Xj. The T a are fundamental representation SU(A^ C ) color matrices 
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normalized so that Tr(T a T b ) = 5 ab . The color-ordered amplitude A^ ee is invariant under a 
cyclic permutation of its arguments. 

We describe the amplitudes using the spinor helicity formalism. In this formalism ampli- 
tudes are expressed in terms of spinor inner-products, 

(jl) = (j-\l + )=u4k J )u + (k l ), [jl] = (j + \r) =«+(*;>_(*,), (2) 

where u±(k) is a massless Weyl spinor with momentum k and plus or minus chirality [16, 18]. 
Our convention is that all legs are outgoing. The notation used here follows the standard 
QCD literature, with 

(ij)\ji] = 2k i -k j = s ij . (3) 

(Note that the square bracket [ij] differs by an overall sign compared to the notation com- 
monly used in twistor-space studies [1].) For non-MHV amplitudes we also use the abbrevi- 
ated notation, 

(i + \(a + b)\j + ) = (i + \{L + h)\j + ) , 
(i-\(a + b)(c + d)\j + ) = + + ■ (4) 

We denote the sums of cyclicly-consecutive external momenta by 

Kl^ EE (h + k i+1 + ■■■ + kj-i + kjY , (5) 

where all indices are mod n for an n-gluon amplitude. The invariant mass of this vector is 
Sj...j = Kf -. In the seven-point case, using momentum conservation we just need to consider 
two- and three-particle invariant masses, which are denoted by 

s ij = [ki + kj) = 2ki ■ kj, Sijk = {k% + kj + kk) ■ (6) 

In color-ordered amplitudes only invariants with cyclicly-consecutive arguments appear, 
and Sj,j+i i j+2- 

The simplest of the partial amplitudes are the maximally helicity-violating (MHV) Parke- 
Taylor tree amplitudes [7] with two negative-helicity gluons and the rest of positive helicity, 



A tree MHV/i o 

A ^ k U ' 2 '---' nj "'(12) (23>---(nl>' 
where j and k label the negative-helicity legs. 



For one- loop amplitudes, the color decomposition is similar [19]. When all internal par- 
ticles transform in the adjoint representation of SU(iV c ), as is the case for M = 4 supersym- 
metric Yang-Mills theory, we have 

|n/2j+l 

A 1 n - loo »({k u \ u a l })= E Gr„ ;c (a) A n;c (o-) , (8) 

where [^J is the largest integer less than or equal to x. The leading color-structure factor 

Gr n;1 (l) = N C Tr(T ai ---T a "), (9) 

is N c times the tree color factor. The subleading color structures are given by 

Gr n;c (l) = Tr(T ai • • -T"- 1 ) Tr(T ac • ■ -T a "). (10) 

S n is the set of all permutations of n objects, and S n;c is the subset leaving Gr n;c invariant. 

The one-loop subleading-color partial amplitudes are given by a sum over permutations 
of the leading-color ones [9], 

A n;c (l,2,..., C -l;c,c + l,...,n) = (— l) c_1 ^ A^a) , (11) 

aeCOP{a}{/3} 

where a { e {a} = {c-l,c-2, . . . , 2, 1}, fa e = {c,c+l, . . . ,n-l,n}, and COP{a}{(3} 
is the set of all permutations of {l,2,...,n} with n held fixed that preserve the cyclic 
ordering of the ctj within {a} and of the within {/?}, while allowing for all possible relative 
orderings of the ctj with respect to the fa. To obtain the full amplitude, therefore, we need 
only compute the leading-color single-trace partial amplitudes. The simple relation (11) 
between leading- and subleading-color contributions is special to one loop; at higher loops, 
new non-planar structures enter the subleading-color partial amplitudes. 

We also collect here the results for the MHV partial amplitudes in the M = 4 theory [9] . 
These amplitudes appear in various kinematic limits of the NMHV amplitudes, so it is 
useful to understand their structure. The MHV amplitudes are a simple linear combination 
of certain box integral functions, 

<r 4MHV = crir xV*. (12) 
The factor V* (n > 5) depends on whether n is odd (n = 2m + 1) or even (n = 2m), 

m— 1 n 

(/i 2 )~ € Vf m+1 = ^ ^ F 2me (Si...( i+r ), S(i-i)...(j +r -i), Si...( i+r ._i), S(i +r+ i)...(i_2)) 
r=2 j=l 
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+ F lm ( s i-3,i-2, Si-2,i-l, Sj...(i-4)) , 

i=l 
m—2 n 

(^ 2 ) ^2™ — ^ ^ ^ 2me ( s j...(j+r), S(j-l)...(i+r-l), Sj...(j +r _l), S( i+r+ i)...(j_2)) 

r=2 i=l 
n 

+ ^ -^ lm (Sj-3,j-2, Sj_2,i-1, Sj...(i_4)) 
i=l 
n/2 

i...(i+m-l)) s (i-l)...(t+m-2)) Si...(j+ro-2) , S(i+m)...(i-2 ))• (13) 

i=i 

The box integral functions F are defined in appendix A. They are essentially scalar box 
integrals, but multiplied by convenient normalization factors in order to remove all power- 
law behavior in the kinematic invariants, leaving only logarithmic behavior. By dimensional 
analysis, the scale \i enters all dimensionally-regulated (unrenormalized) one-loop amplitudes 
as an overall factor of (/x 2 ) e . 

The four-point and five-point amplitudes appear in the multi-particle factorization limits 
of the seven-point amplitudes, as discussed in section V and appendix B. The n = 4 case is 
given by 

<f 4 = 2 c r (u 2 ) £ AT e x F 0m ( Sl2 , s 23 ) , (14) 
for all non-vanishing helicity choices. For n — 5, the expression (13) reduces to 

<f 4 = CP (/! 2 ) £ Af C [F lm ( Sl2 , S 23 , S 45 ) + F lm (.S 23 , S 3 4, S 5 l) + F lm ( S3 4, S 45 , S 12 ) 

+F lm (s 45 , s 5 i, s 23 ) + F lm (s 5 i, s 12 , s 34 )] • (15) 

Besides the Af = 4 MHV amplitudes, a number of other infinite series of one-loop am- 
plitude have been computed. The n-point MHV gluon amplitudes with an M = 1 chiral 
multiplet in the loop were also computed using the unitarity method. In QCD one-loop 
n-point amplitudes with identical helicities are also known [42]. At two loops less is known. 
A conjecture for the planar two-loop MHV amplitudes in terms of one-loop amplitudes was 
presented in ref. [12], suggesting that the simplicity uncovered at one loop persists to higher 
loops in the planar, leading-color limit N c — > oo. 

The one-loop non-MHV six-point amplitudes in M — 4 super- Yang-Mills theory were 
computed in ref. [10]. Here we will compute all non-MHV seven-point amplitudes. (For 
n = 6 and n = 7, 'non-MHV coincides with 'next-to-MHV'.) We expect that these new am- 
plitudes will be helpful for unraveling the full n-point twistor-space structure. Indeed, their 

10 



structure provided important guidance for obtaining the coefficients of a class of three-mass 
box integrals appearing in the all-n next-to-MHV amplitudes, which we present in section VI. 

III. CALCULATIONAL APPROACH 

It seems clear that calculating the seven-point amplitude by computing the 227,585 
contributing Feynman diagrams one by one is not the best way to proceed. And an all- 
multiplicity result is simply not accessible to this standard approach. 

Fortunately, the unitarity-based technique is better for this application. (For a recent 
review, see ref. [26].) The basic idea is to reconstruct a one-loop amplitude from its cuts 
or absorptive pieces, which are products of tree amplitudes. In general, one must calculate 
these cuts in D dimensions. In this case there are no ambiguities in reconstructing ampli- 
tudes in any massless theory. In the special case of supersymmetric theories, one can even 
evaluate the cuts 'in four dimensions' — that is, by assigning four-dimensional helicities to 
the states crossing the cut — without encountering any ambiguities, for all terms in the am- 
plitude which survive as e — > [10]. This property reflects the better ultraviolet behavior of 
these theories. Then the one-loop cuts reduce to products of tree-level helicity amplitudes. 
Starting from tree amplitudes rather than diagrams means that the extensive cancellations 
that occur in gauge theories are taken into account before any loop integrations are done, 
which greatly reduces the complexity of the calculations. Because infinite series of tree- 
level amplitudes are known, it also opens the door to calculating infinite series of one-loop 
amplitudes. 

The reconstruction is done by identifying the corresponding Feynman integrals whose cuts 
yield the original integrands. This is easier than doing the dispersion integrals explicitly, 
because we take advantage of an extensive machinery for simplifying and computing such 
integrals. 

In calculating an n-point amplitude in gauge theory, one expects to encounter n-point 
tensor integrals with tensor rank up to n. In supersymmetric theories, the maximal tensor 
rank is lower, a reflection of the better ultraviolet behavior of the theory [10]. In M = 4 
theories, the maximal tensor rank is reduced to n — 4 [9, 43]. In addition, thorough use of 
the spinor-helicity representation of external gluon polarization vectors can reduce both the 
tensor rank and multiplicity of integrals we encounter upon reconstructing them from the 
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cuts. 

Nonetheless, we do have to treat some higher-point and tensor integrals. At one loop, on 
general grounds, one can reduce any loop integral in D dimensions to a combination of box 
integrals, triangles, and bubbles. In Af = 4 supersymmetric theories, only scalar box inte- 
grals are needed, and no triangles or bubbles. The tensor higher-point integrals encountered, 
where loop momenta appear in the numerator of the integrand, could be reduced to scalar 
ones via brute- force Brown-Feynman/Passarino-Veltman techniques [44], or alternatively 
via Feynman-parameter methods. However, we have found it useful to follow a method 
adapted from that originally proposed by van Neerven and Vermaseren [45], relying on the 
fact that any vector in four dimensions can be written in terms of a basis of four vectors. 
For example, we may expand the vector as 

where e' ip2PSP4 is shorthand for & ^apPzPzPii an d is the Levi-Civita tensor. Of course, 
we are working with dimensionally-regulated integrals, but we can split the loop momentum 
vector £ M into two parts, the four-dimensional components and the (— 2e)-dimensional ones. 
The four-dimensional components can be expanded according to eq. (16), where the pi will 
be chosen to be either external momenta or sums of external momenta (see below). In this 
case, we can take advantage of the fact that 

= l(t 2 -(t-P*) 2 +Pl), (17) 

in order to cancel propagators, or reduce the tensor rank of the integral. In the unitarity- 
based method, £ 2 is taken to be zero, and so when we are using the spinor-helicity represen- 
tation, what appears in the numerator are only dot products of the loop momentum, and 
spinor 'sandwiches' of the loop momentum between spinors representing external momenta. 
These objects depend only on the four-dimensional components of the loop momentum. (Be- 
cause we are evaluating the cuts in four dimensions, we have dropped some terms in the cut 
integrand which are proportional to the square of the (— 2e)-dimensional components. These 
pieces lead only to (9(e)-suppressed terms in the amplitudes, for supersymmetric theories.) 

We do not choose a fixed basis of external vectors for all expansions. For terms containing 
multiple factors of the loop momentum in the numerator, we proceed in stages, expanding 
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the successive occurrences of the loop momentum one by one. At each stage, we choose the 
basis set for the expansion to consist of external momenta, not of the amplitude, but of the 
subject integral. In particular, after the first stage of this procedure, one of the external 
legs will consist of a sum of the original external momenta, because of the cancellation of 
one of the propagators using eq. (17). We choose such a sum to be an element of the basis 
set for the second stage. Using the new basis set in the second stage ensures that the loop- 
momentum-containing quantities £ 2 and {£—pi) 2 appearing on the right-hand side of eq. (17) 
will continue to cancel some propagator in the integral at this stage. We proceed similarly 
in the further stages. We call this approach a 'pivoting' technique. Since dot products of 
the loop momentum with external legs of any one-loop integral are always expressible in 
terms of inverse propagators and external invariants, as in eq. (17), our choice always allows 
a given integral to be reduced into a sum of lower-rank and lower-point integrals. (Were we 
to have chosen a fixed external basis set, expansions beyond the first would not necessarily 
be expressible in terms of inverse propagators of the daughter integrals, complicating the 
analysis.) 

In principle the pivoting technique can be applied to any cut (or even to full uncut 
diagrams if desired). For the Af = 4 application, however, it is only necessary to consider 
cuts in three-particle channels, those with invariants crossing the cut. That is 

because box integrals contain cuts in several channels, and for the seven-point process at 
least one channel is always a three-particle channel. Such a cut is the product of a five-point 
tree amplitude and a six-point tree amplitude, or it may be a sum of such terms, the sum 
corresponding to multiple members of the M = 4 supermultiplet crossing the cut. It is 
useful to divide the cuts into 'singlet' and 'nonsinglet' pieces. The singlet pieces have a total 
helicity of ±2 crossing the cut. They can only come from gluonic intermediate states. The 
nonsinglet pieces have a total helicity of crossing the cut, and receive contributions from 
the entire M = 4 supermultiplet, gluons, fermions and scalars. One can further divide the 
singlet and nonsinglet pieces according to whether the six-point tree amplitude(s) are MHV 
or non-MHV. So there are four types of cut building blocks: 

• singlet-MHV. This cut is the simplest type to evaluate, and can be done directly in a 
compact form for all n; see ref. [40] and section VI). 

• nonsinglet-MHV. This cut is also simple; the sum over Af = 4 states can be performed 
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FIG. 1: A generalized 'triple cut'. The three propagators cut by the dashed lines are required to 
be 'open'. 

as in ref. [10], resulting in an expression which is identical to the singlet-MHV case, 
except for overall, loop- momentum-independent, prefactors. 

• singlet-non-MHV. This cut is more complicated, but still only the product of two 
helicity amplitudes. 

• nonsinglet-non-MHV. This cut is the most complicated, but only needed to be evalu- 
ated as a cross-check. 

It is also possible to consider generalized cuts which provide very useful additional infor- 
mation about loop amplitudes [26, 46]. The example of a 'triple cut' is shown in fig. 1. It is 
not a cut in the traditional sense of having a phase-space integral which generates the imagi- 
nary part in a single channel. It corresponds to a collection of underlying Feynman diagrams 
which have three propagators 'open', or uncancelled, in the same way that an ordinary cut 
requires two propagators to be open. Terms which vanish as the 'cut' propagators go on 
shell may be ignored. The contributing terms can be represented as the product of three 
tree amplitudes. We will see in section VI that analysis of such cuts provides additional 
information, often quite easily. 

Application of the pivoting technique to each type of cut yields an expression in terms 
of external spinor invariants, and scalar integrals up to n points. To reduce the six- and 
seven-point scalar integrals we used the formulae from appendix VI of ref. [47]. To reduce 
the five-point scalar integrals we made use of two different reduction formulae. The first 
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formula is also from ref. [47] in the equivalent form from ref. [48], 

n 

^ = ^E c ^' ) + ec o^ 6 - 2e , (18) 

where 1^ is the daughter box integral derived from the pentagon by omitting (the negative 
of) the propagator between legs (j — 1) and j; and where 

n 

1=1 

n 

co = $>, (20) 
i=i 

in terms of SV,- = Sji = —Kf (j_i)/2 (and Su = 0). The pentagon integral in 6 — 2e 
dimensions, I^ =6 ~ 2e , is finite as e — > 0, and it has a manifest prefactor of e in eq. (18), so 
its contribution may be neglected to the order we are working. 

These reductions are valid in dimensional regularization. The six- and seven-point reduc- 
tion formulae are equivalent to those of Melrose [49] , and van Neerven and Vermaseren [45] . 
The five-point reduction formula is equivalent after dropping the 0(e) terms. The five-point 
reduction formula (18) is valid for integrals, but does not hold point-by-point for the in- 
tegrands. We found it convenient for certain parts of the calculation to use a reduction 
formula which is valid point-by-point. We can obtain such a formula by adding appropriate 
terms containing the Levi-Civita tensor and linear in the loop momentum, e M ' £ M . . ., which 
vanish upon performing the loop integration. The advantage of including these terms in the 
reduction is that one can check all algebra numerically, directly on the integrand. 

Once we have reduced the integrand to a sum of four-point integrands, we are done, 
because amplitudes in the maximally supersymmetric theory can be expressed in terms of 
scalar boxes, with no triangle or bubble integrals. (Due to Af = 4 supersymmetry, each term 
in the cut integrand with r — 2 non-cut propagators starts out with no more than degree 
r — 4 tensor integrals. Given this starting point, the pivoting reduction ensures that tensor 
boxes never appear.) 

This procedure gives us the coefficients of the standard box integral functions. The forms 
emerging directly from the reductions are generally quite complicated. However, we can 
simplify them further, as discussed below. Somewhat surprisingly, many of them have a 
remarkably simple form. Indeed, the most 'complicated' integrals that appear, the three- 
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mass boxes, have single-term coefficients, and the most complicated coefficients have only 
four terms. 

The coefficients have various momentum-space singularities. Some of the coefficient sin- 
gularities are also singularities of the full one-loop amplitude, such as the collinear singu- 
larities or poles in multi-particle invariants s^-fc. Other singularities are spurious. These are 
introduced by the reduction procedure. Some of these singularities disappear coefficient by 
coefficient; others survive. They will cancel in the amplitude as a whole, as we will discuss in 
section V; but this cancellation involves the box integrals and their behavior in an essential 
way. That is, each coefficient can (and does in general) possess singularities not present in 
the amplitude as a whole. The pentagon reduction (18), in particular, gives rise to denom- 
inators which induce such singularities. They can be read off from eq. (19) by inspecting 
the determinant of the appropriate pentagon SV,- matrix. In the seven-point calculation, we 
encounter one-mass and two-mass pentagons. The former can be treated as a special case of 
the latter. The non-adjacent two-mass pentagon (with external legs (ki, k 2 , P, k 3 , Q), capital 
letters denoting massive legs) gives rise to the following denominator, 

2k 1 -k 2 ((P + k 2 ) 2 (P + k 3 ) 2 -P 2 (P + k 2 + k 3 ) 2 )((Q + k 3 ) 2 (Q + k 1 ) 2 -Q 2 (Q + k 1 + k 3 ) 2 ). (21) 

Each of the last two factors can be further factorized using the spinor identity [46], 

(P + h) 2 (P + k 2 ) 2 - P 2 (P + k 1 + k 2 ) 2 = (l-|^|2-)<2-|/>|l-> . (22) 

This type of denominator thus develops singularities when P lies in the (space-time) plane 
spanned by k\ and k 2 , 

P" = ai k1 + a 2 A# , (23) 

for arbitrary values of a\ and a 2 . For this reason, we call the factor in eq. (22) and the 
corresponding singularity a 'planar' or 'back-to-back' one. 

The adjacent-mass two-mass pentagon (with external legs (ki, k 2 , k 3 , P, Q)) gives rise to 
the denominator, 



Ak v k 2 k 2 -k 3 



{k 1 +k 2 +k 3 ) 2 {Q+k 1 ) 2 {P+k 3 ) 2 -2k v k 2 P 2 {Q+k 1 ) 2 -2k 2 -k 3 Q 2 {P+k 3 ) 2 j. (24) 

The last factor, in brackets, also has a spinor factorization (for k\ + k 2 + k 3 + P + Q = 0), 

(l-|^|3+><3 + |^|l->. (25) 
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It therefore vanishes whenever 



P = aiki + a 2 k 2 + a 3 k 3 , 



(26) 



with 



ki ■ k 2 + (1 + a 3 )ki ■ k 3 



(27) 



a 2 = — ai 



aifci • fc 2 + (1 + a 3 )/c 2 • h ' 



To distinguish it from the other kind of denominator and singularity, we will call this type 
'cubic'. 

In computing the all-n coefficient discussed in section VI, one also encounters denomina- 
tors arising from the reduction of certain three-mass pentagon integrals. It turns out that 
they also can be expressed in terms of these sorts of 'cubic' factors. 

Knowledge of all of the singularities of the box coefficients allows one to simplify them 
'numerically'. The basic idea is to generate kinematic points numerically which are close to 
the various singular regions: nearly collinear kinematics for each of the seven color- adjacent 
pairs (i, i + 1); kinematic points where each of the seven multi-particle invariants 5^+1,1+2 
almost vanish; and points near the relevant 'planar' and 'cubic' singularities. By studying 
the behavior of the large analytical expressions for the box coefficients in these regions, one 
can write down all the allowed factors in the denominator of the expression. In the collinear 
limit where ki is nearly parallel to /cj+i, it is helpful to be able to rotate ki and about 
their common sum. The complex phase behavior under this rotation reveals whether the 
denominator contains (i (i + 1)) or [i (i + 1)] (or both). 

The numerator of the box coefficient can often be completely determined by the spinor 
homogeneity properties, which dictate how many net powers of vs. occur, in terms 
of the helicity of leg i. In a few cases, we solved for a coefficient by writing down a sum 
of all possible expressions of the right dimension and spinor homogeneity, and comparing it 
numerically to the actual expression at a number of randomly-generated phase-space points. 

The easy-two-mass box coefficients are the simplest to determine in this way, because 
they cannot have any 'cubic' singularities, and because they turn out to be composed of 
just one term. (The easy-two-mass box integral has a massless leg interposed between each 
massive one, which prevents it from being obtained from an adjacent-mass pentagon by 
cancelling a propagator. The adjacent-mass pentagons are the only source of cubic singu- 
larities.) For example, in examining the coefficient C147 for A^= 4 (l~, 2~, 3~, 4 + , 5 + , 6 + , 7 + ), 
one quickly finds that it vanishes strongly as legs 2 and 3 become collinear, and as legs 
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6 and 7 become collinear. The factors of (2 3) 3 [6 7] 3 in the numerator of eq. (41) can be 



established in this way. The denominator factors in eq. (41) are all easily obtained numer- 
ically as well, except for the question of whether (1 + | (6 + 7) |5 + ) or its complex conjugate 
(1~| (6 + 7) |5~) correctly describes the 'planar' singularity, and similarly for (6 + | (7 + 1) |2 + ) 
vs. (6~| (7 + 1) |2~). However, only one of these four possibilities can give the right spinor 
homogeneity for legs 1, 2, 5 and 6. (The possibility of additional factors in the numerator 
is excluded in this case by dimensional analysis.) 

In addition to simplifying the coefficients obtained from a direct calculation, we can also 
employ another tool. The structure of infrared singularities provides equations which can 
be used as consistency checks, as we shall discuss in section V, or alternatively to solve for 
some of the coefficients. The infrared singularities in the amplitude are known on general 
grounds [50] to be, 



That is, the 1/e poles are proportional to the seven-gluon NMHV tree amplitude [51] and 
contain only logarithms of nearest-neighbor two-particle invariants. This implies that the 
coefficient of any given ln(— s iji+1 )/e must be equal to the tree; and the coefficient of any 
ln(-s i,i+i,j+2)/e must vanish. On the other hand, the scalar box functions whose coeffi- 
cients we are computing contain both of these sorts of terms, and so both types of equa- 
tions are non-trivial. Each box function in eqs. (A2)-(A5) contains various ln(— s^+^/e 
and ln(— Sj ii+li j +2 )/e terms with coefficients 0, ±1 or ±|. Thus the constraints arising 
from eq. (28) become simple linear relations among the coefficients, some of which involve 
the tree amplitude. (Other linear relations may be determined from numerical evaluation of 
the box coefficients.) We use the 1/e pole information as part of the simplification process, 
but the simplest forms of the coefficients do not satisfy it manifestly. That is, we shall obtain 
alternate, compact representations for the tree amplitudes, by re-imposing the constraints 
from eq. (28) after completing the simplifications. 

It is also possible to relate box coefficients from one of the four seven-point NMHV helicity 
amplitudes to those from the other three, using properties of the tree amplitudes which enter 
the cuts. Such relations are possible any time a cut that is sensitive to a given box integral 
does not contain the 'nonsinglet-non-MHV type of contribution. Under such circumstances, 
there are supersymmetry Ward identities (SW) [21] which can be used to permute around 




(28) 
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S-NMHV 123 = Af c ((-h) + , 1+, 2", 3", l\) A*° c ((-£ 4 )", 4" 5+, 6+, 7+ £7) . (30) 



the helicities on one side of the cut. 

Consider for example the S123 cut of A^= 4 (l~, 2~, 3+, 4~, 5 + , 6 + , 7 + ). It has a singlet- 
non-MHV term, which is given by 

S-NMHV 123 = i A'rd-iiV, I", 2-, 3+, 4) Af c ((-£ 4 )-, 4", 5 + , 6 + , 7+, £7) , (29) 
integrated over phase space. Using a SWI on the left-side of the cut, the cut becomes 

"<2 3)< 

Except for the spinor-product prefactor, this expression is equivalent to the singlet-non- 
MHV piece of the s 234 cut of A^= 4 (l~, 2", 3", 4+, 5+, 6+, 7 + ), after the relabelling 1 <-> 4, 
2 <-> 3, 5 <-> 7 (and an overall minus sign). On the other hand, the nonsinglet-MHV term in 
the above S123 cut is given by 

NS-MHV 123 = % ^ 5 ree ((-4)- A , I", 2", 3+, £ A ) ^ ec ((-^)" A , 4-, 5+, 6+, 7 + , £ A ) 

A 

(3+1(1 + 2)14+) t A + A , 

" V|(2 + 3)|4+) 4 V 5 (( } ^ 

xAr(M4)- A ,4-,5 + ,6+,7+,^). (31) 

Again, except for the spinor-product prefactor, this expression is equivalent to the nonsinglet- 
MHV piece of the s 2 34 cut of A^ 4 (l~, 2", 3", 4+, 5+, 6+, 7+), after applying the same rela- 
belling and minus sign. 



IV. RESULTS 

In this section we present the results for the four independent one-loop 7-point NMHV am- 
plitudes in M = 4 SYM: A^=\l~, 2" 3", 4+ 5+ 6+, 7+), A^=\\- , 2" 3+ 4", 5+, 6+, 7+), 
A^' 1 = 4 (l-,2-,3+,4+,5-,6+,7+), and A^f 4 (l-, 2+, 3", 4+, 5", 6+, 7+). All other A^ A 
NMHV helicity amplitudes can be obtained from these by the three operations: 

1. parity, which exchanges + and — helicities and is implemented on the basic spinor 
products by the action (ij) <-> [j i], 

2. reflection symmetry, essentially charge conjugation invariance, which states that 
A^r 4 (l Al ,2 A V..X™) = (-l) n ^i =4 (™ A V--,2 A2 ,l Al ), and 
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3. cyclic symmetry, since each A^ 4 is the coefficient of a cyclicly-invariant color trace. 

The subleading-color partial amplitudes appearing in eq. (8), A^ c =4 for c > 1, are obtained 
by sums over permutations of the leading-color A^f 4 [9], using eq. (11). 

The M = 4 SYM amplitude is expressible as a sum of scalar box integrals I4, or equiva- 
lently the box functions F given in appendix A. For the color-ordered amplitudes Ay ; i, the 
kinematics of each box integral can be obtained from the heptagon diagram with external 
legs in the order 1,2,3,4,5,6,7, by deleting three different propagators, say k. We label 
each box function F, and each kinematic coefficient c multiplying it, by this triplet of in- 
tegers. However, to avoid confusion with a twistor-space 'line operator' F^ to appear in 
section VII, we call these box functions B(i,j,k) instead of F(i,j,k). Thus we write the 
leading-color partial amplitude as 

A«= 4 = tc r (fi 2 Y ]T c ijk B(i,j, k) , (32) 

where 

1 r(i + e)r 2 (i- e ) 

Cr " (4^ f(T—2e) (33) 
is a ubiquitous prefactor. We label the propagators so that in the heptagon diagram the 
external leg labeled p lies between propagators labeled p and p + 1 (mod 7). Thus in the box 
B(i,j, k), leg % is grouped into a single massive leg with leg i — and similarly for legs j and 
k. If two members of the set {i,j,k} are cyclicly adjacent, the corresponding massive leg 
of the box integral contains three massless legs. If all three of {i,j, k} are cyclicly adjacent, 
that massive leg contains four massless legs. Some examples of this labeling are shown in 
fig. 2. 

The kinematic arguments of the boxes are s and t, the squares of the sums of momenta 
emerging from two adjacent vertices of the box, followed by the 'masses', the squares of the 
sums of momenta emerging from individual vertices of the box (when these are nonzero). 
The easy two-mass box has two reflection symmetries, so it doesn't matter which invariant 
is s and which t. For the hard two-mass box, s is the invariant formed by the two adjacent 
massless legs. For the three-mass box, s is the invariant formed by the massless leg and the 
first massive leg following it cyclicly. 

For a given helicity amplitude, the number of box functions, and box coefficients, is the 
number of un-ordered integer triplets (i,j,k), where each integer runs from 1 to 7, and all 
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(a) (b) (c) (d) 

FIG. 2: Examples of box integral functions B(i,j,k) appearing in 7-point amplitudes; the ar- 
guments k are circled: (a) the one-mass box 5(1,2,7) = i ?lm (s34, S45, S345), (b) the 'easy' 
two-mass box 5(1,2,5) = F 2rae (s^, S456, S45, S712), (c) the 'hard' two-mass box 5(1,2,4) = 
F 2mh (s 56 , 5345,^712,534), and (d) the three-mass box 5(1,3,5) = F 3m (s67i, S456, 571, S23, £45)- 

three are unequal. This number is just (3), or 35. Another way to count this number is to 
observe that there are 5 different kinds of box in the 7-point case: 

• the one- mass box shown in fig. 2a, plus cyclic permutations (7 boxes), 

• the easy-two-mass box shown in fig. 2b, plus cyclic permutations (7 boxes), 

• the hard-two-mass box shown in fig. 2c, plus its 'reflection' in which the 2-leg cluster 
precedes the 3-leg cluster, plus cyclic permutations (14 boxes), 

• the three-mass box shown in fig. 2d, plus cyclic permutations (7 boxes). 

In the remainder of this section we present the coefficients c^fc of all the box functions, 
using symmetries whenever possible. In section VII we shall describe the twistor-space 
structure and other general properties of these coefficients, which fit nicely into a uniform 
pattern. 

A. A^= 4 (l-,2-,3-,4+,5+,6+,7+) 

This partial amplitude possesses a reflection, or "flip" symmetry. We define the flip 
operation (for this partial amplitude) by 

X| flip = -X(l <-> 3,4 <-> 7,5 <-> 6) , (34) 
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i.e. a two-fold exchange of labels accompanied by an overall minus sign. The tree amplitude 
obeys 

4 rcc (l-, 2" 3-, 4+, 5 + , 6+ 7+) = 4 rcc (l-, 2', 3", 4+ 5+ 6+, 7+) , (35) 

flip 

as does the one-loop amplitude, 

<f 4 (1~, 2", 3", 4 + , 5 + , 6 + , 7+) = <r 4 (l", 2", 3", 4+ 5+, 6+, 7+) . (36) 

flip 

We can write all of the box coefficients in terms of the following quantities: 

-(2 3) 3 (5+|(6 + 7)|l+) 3 

A ~ S2 34 S 67i(34)(67)(71)(5 + |(3 + 4)|2+)(6-|(7+l)(2 + 3)|4 + ) ' 1 ) 

r -(3-|(4 + 5)(6 + 7)|l+) 3 

B ^345^671 (3 4) (4 5) (6 7) (71) (2+| (3 + 4) |5+) (2+| (7+1) |6+> ' W 

(2 3) 3 (7+|(5 + 6)|4+) 3 

Cl36 ~ (3 4) (4 5) (5 6) [71] (1+| (2 + 3) |4+> (6"| (7 + 1)(2 + 3) |4+> (2"| (7 + 1)(5 + 6) |4+) ' 

(39) 

(12) 3 



0146 [56] X (34) (5+| (3 + 4) |2+> (2-| (3 + 4) (5 + 6) |7+> 

< 23 > 3 (40) 



(71)(6+|(7+l)|2+)(2-|(7+l)(5 + 6)|4+) ' 

r = (2 3) 3 [6 7] 3 

147 " sen (3 4) <4 5> [71] <1+| (6 + 7) |5+> (6+ 1 (7+ 1) |2+> ' 1 ) 

(5-|(6 + 7)(l + 2) |3+) 3 

247 " [12](34)(45)(56)(67)(l + |(6 + 7)|5 + )(2 + |(3 + 4)|5 + )(7-|(l + 2)(3 + 4)|5 + ) ' 

(42) 

„ = (2 3) sl 67 

346 " S2 34 (34)(56)(67)(l + |(2 + 3)|4 + )(l + |(6 + 7)|5 + )(2-|(3 + 4)(5 + 6)|7 + ) ' 1 } 
r = (4+|(2 + 3)|l+) 3 

347 ' S23 4 [23][34](56)(67)(71)(2+|(3 + 4)|5 + ) ' 1 1 

° 567 = [1 2] [2 3] (4 5) (5 6) (6 7) (1+| (2 + 3) |4+) (3+| (1 + 2) |7+) ' (45) 

Note that cb and are flip symmetric. 

The quantities (37)-(45), plus their images under the flip operation, are not all indepen- 
dent. They satisfy 

[C136] [flip + C 3 46 + C 5 67 = C 24 7 + C 3 47 , (46) 
Ci4 6 + C M7 + C 346 = C A + Ci 36 , (47) 
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and 



Ay CC = [C A + Ci36]|fiip + Ci47 + C346 + C567 (48) 
= C B + C 347 + [C 3 4 7 ]| flip . (49) 

Equation (49) in particular provides a compact representation of the tree amplitude which 
has the proper collinear behavior manifest in all channels, resembling some "twistor-space" 
constructions [23]. Equation (49) is manifestly flip symmetric, whereas eq. (48) is not. 
In terms of these quantities, the remaining box coefficients are given by 



and 



C235 


— c 237 — c 256 — c 257 = c 357 = c 367 = 0, 


(50) 


C236 


= C 2 67 = C356 = C 56 7, 


(51) 


Cl34 


= c A + cu 7 , 


(52) 


cm 


= c A + c 347 , 


(53) 


Cl67 


= C A + Ci 36 + C347, 


(54) 


C234 


= C247 + C347, 


(55) 


C345 


= C B + C U7 + C347, 


(56) 


C457 


= C B + C147, 


(57) 


C467 


= C346 + C347, 


(58) 




c 123 — C234 1 flip , 


(59) 




C124 = C134 flip , 


(60) 




Cl25 — c 347 flip j 


(61) 




Cl26 — C346 flip , 


(62) 




Cl27 — C345 flip , 


(63) 




c 135 — c 247 flip j 


(64) 




C145 = C147 flip , 


(65) 




Cl56 — Q67 flip , 


(66) 




Cl57 = C457 flip , 


(67) 




C245 — C137 flip , 


(68) 
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C246 — c 136 Iflip 5 

(69) 

C456 = c 167 | flip . (70) 
The flip relations can also be summarized as 

°ijk = %fc|flip, (71) 

where 

i = (5 - i) mod 7, (72) 
and ijk should be written in ascending order. 



B. 4r 4 (l-,2-,3+,4-,5+,6+,7+) 

This partial amplitude is alone among the NMHV seven-point amplitudes in having no 
flip symmetry. The box coefficients are given in terms of the following quantities: 

(24) 4 (5+|(6 + 7)|l+) 3 

A " S671 s 2 34 (23)(34)(67)(71)(5+|(3 + 4)|2+)(6-|(7+l)(2 + 3)|4 + ) ' 1 ' 

r (12) 3 (7+|(5 + 6) |4+) 3 

B Sl2 3 S 456 (23)(45)(56)(7 + |(l + 2)|3+)(6-|(4 + 5)(2 + 3)|l + ) ' 1 ] 

r = -(4-|(3 + 5)(6 + 7)|l+) 4 

° ^345*671 (3 4) (4 5) (6 7) (7 1) (2+ 1 (3 + 4) |5+> (2+ 1 (7 + 1) |6+> (3" | (4 + 5) (6 + 7) 1 1+) ' 

(75) 

(12) 3 (6+|(3 + 5) |4+) 4 

£345*712 (34) (45) (7 1) (6+| (4 + 5) |3+) (6+| (7 + 1) |2+) (5"| (3 + 4)(1 + 2) |7 + ) ' 

(76) 

(12) 3 (3+|(5 + 6) |4+) 4 

^456*712 (4 5) (5 6) (71) (3+| (4 + 5) |6+) (3+| (1 + 2) |7+) (2"| (7 + 1)(5 + 6) |4+) ' 

(77) 

(7 + |(l + 2)|4+) 4 
s 712 [1 2] (3 4) (4 5) (5 6) [7 1] (2+ 1 (7 + 1) |6+) (7+ 1 (1 + 2) |3+) 
((3 + |(6 + 7)|l + )(46) + (3+|5|4 + )(16)) 4 



CD = 



CE 



C 125 - - n 01 10 a\ 1 a T\ r-711 /0+1 m , 7\ h , io+\ > 



° 135 ~ [23](45)(56)(67)(71)(3 + |(4 + 5)|6 + )(2 + |(7+l)|6 + ) 

X (6-| (7 + 1)(2 + 3) |4+) (6- 1 (4 + 5) (2 + 3) |1 + ) ' (79) 
(2 4) 4 (7+|(5 + 6)|4+) 3 

C136 " 



(23) (34) (45) (56) [7 1] (1+| (2 + 3) |4+) 



X (2-|(7+l)(5 + 6)|4+)(6-|(7+l)(2 + 3)|4 + ) ' (80) 
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(1 2) 3 [3 5] 4 

° 145 = S345 [3 4] [4 5] (6 7) (7 1) (3+| (4 + 5) |6+> (5+| (3 + 4) |2+> ' (81) 



(1 2) 3 (2 4) 4 [5 6] 3 

Cl46 — 



(2 3) (3 4) (7 1) (5+| (3 + 4) |2+> (6+| (7 + 1) |2 + ) 

1 

X (2-| (3 + 4) (5 + 6) |7+) (2- 1 (7 + 1)(5 + 6) |4+) ' (82) 

c = (2 4) 4 [6 7f 

0147 S6n (23)(34)(45)[71](l+|(6 + 7)|5+)(6+|(7+l)|2+) ' ^ 

C236 = (3 + l(l + 2)|4+) 4 (84) 

236 Sl23 [12][23](45)(56)(67)(l + |(2 + 3)|4+)(3 + |(l + 2)|7 + ) ' 1 1 

r = ( 12 ) S 567 /nr\ 

237 " Si 2 3 (23)(56)(67)(4+|(2 + 3)|l + )(4+|(5 + 6)|7+)(3-|(l + 2)(6 + 7)|5+) ' K ™> 

(12) 3 (3+|(5 + 6)|7 + ) 4 

[34] (56) (6 7)(71)(3+|(l + 2)|7+)(4+|(5 + 6)|7+) 

X (2"| (3 + 4) (5 + 6) |7+) (5- 1 (3 + 4)(1 + 2) |7+) ' (86) 
(4-|(l + 2)(6 + 7)|5+) 4 

C247 ~ 



C246 



[1 2] (34) (4 5) (5 6) (6 7) (1+| (6 + 7) |5+) (2+| (3 + 4) |5+) 



X (7-|(l + 2)(3 + 4)|5+)(3-|(l + 2)(6 + 7)|5+) ' (8?) 



(12) 3 [5 6] 3 

c 256 



C257 = 



S4 56 (2 3) [45] (7 1) (4+| (5 + 6) |7+) (6+| (4 + 5) |3+) ' 

(12) 3 (34) 3 [6 7] 3 

(23) (45) (7+|(l + 2)|3+) (6+| (4 + 5) |3+) 

X (5-| (6 + 7)(1 + 2) |3+) (1-| (6 + 7) (4 + 5) |3+) ' (89) 

(2 4) 4 S j 67 

S234 (2 3) (34) (5 6) (6 7) (1+| (2 + 3) |4+) (1+| (6 + 7) |5+) (2"| (3 + 4) (5 + 6) |7+) ' 

(90) 

c = (3 + l(2 + 4)|l+) 4 

S234 [23][34](56)(67)(71)(2+|(3 + 4)|5+)(4+|(2 + 3)|l+) ' 1 ) 

(12) 3 (5+| (6 + 7)|l+) 

C357 " 



C346 — 



3 



(23) [45] (6 7) (71) (4+|(2 + 3)|l+) 

1 

X (3-|(4 + 5)(6 + 7)|l+) (6-|(4 + 5)(2 + 3)|l+) ' (92) 
Again the quantities (73)-(92) are not independent. They satisfy 

CA + C136 = Cue + C147 + C 346 , (93) 
CC + Cl25 + C 2 37 + C 2 57 + C347 = C D + C M7 + C 236 + C 2 46 + C 3 46, (94) 
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CD + C 2 46 + C 2 56 = C E + C145 + C 146 , (95) 

C236 + C 2 46 + C 3 46 = C 237 + C 24 7 + C 3 47, (96) 

CB + C 357 = C 237 + C 256 + C 257 , (97) 

CA + Ci 35 + C M5 = C C + C 34 7 + C 35 7- (98) 

The tree amplitude may be written in terms of the box coefficients in several different ways. 
The shortest expression we have found is: 

A 7 rcc = C B + C C + Ci 25 + c 347 + c 357 . (99) 
The remaining box coefficients are given by 



c l 23 




CE + Ci 3 6 + C 23 6 + C 2 56, 


(100) 


Cl 2 4 




Co + C145, 


(I0l) 


Cl 2 6 




CE + C 25 6, 


(102) 


Cl 2 7 




CC + Ci 25 + C145 + c 257 , 


(103) 


Cl 3 4 




CA + C147, 


(104) 


Cl37 




CA + C 34 7, 


(105) 


c 156 




CE + C125, 


(lUo ) 


Cl57 




CC + C145, 


(107) 


Cl67 




CA + Ci 36 + C 347 , 


(108) 


C 23 4 




C 23 6 + C 2 46 + C 3 46, 


(109) 


C 23 5 




CB + C 25 6, 


(no) 


C 2 45 




CD + Cl25, 


(HI) 


C 2 67 




C 23 6 + c 2 3 7) 


(112) 


C 3 45 




CC + C147 + C 34 7 + C 35 7, 


(H3) 


C 3 56 




CB + c 236 , 


(H4) 


C 3 67 




0, 


(H5) 


Q56 




CD + Ci 2 5 + C 246 + C 256 , 


(116) 


C457 




CC + Cu 7 , 


(H7) 


C467 




C 3 46 + C 3 47, 


(118) 


C567 




CB + c 236 + c 357 . 


(H9) 
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C. ^ 1 = 4 (l-,2-,3+,4+,5-,6+,7+) 

For this partial amplitude we define the "flip" operation, 

X| flip = -X(l <-> 2,3 <-> 7,4 <-> 6), (120) 
which is a symmetry of the tree amplitude, 

4 ree (l-, 2", 3 + , 4+ 5-, 6 + , 7+) = A\ rcc {l-, 2", 3 + , 4+ 5", 6 + , 7+) , (121) 

flip 

and of the one-loop amplitude. The box coefficients c^fc are given in terms of the following 
quantities: 

(12) 3 (7+|(4 + 6) |5+) 4 

° A ~ Sl23S 456 (23)(45)(56)(7 + |(5 + 6)|4+)(7 + |(l + 2)|3 + )(6-|(4 + 5)(2 + 3)|l + ) ' 

(122) 

(12) 3 (4+|(6 + 7) |5+) 4 

S123S567 (2 3) (5 6) (6 7) (4+| (5 + 6) |7 + ) (4+| (2 + 3) |l+> (3"| (1 + 2)(6 + 7) |5 + ) ' 

(123) 

(5-|(3 + 4)(6 + 7)|l+) 4 

5345*671 (3 4) (4 5) (6 7) (71) (2+| (3 + 4) |5+) (2+| (7+1) |6+) 
1 



CB 



cc 



X 



CD 



(3-|(4 + 5)(6 + 7)|l+)' 

(l-|(6 + 7)(3 + 4)|2+) 4 



(124) 



S234S6H (2 3) (3 4) (6 7) (7 1) (5 + | (3 + 4) |2+) (5+| (6 + 7) |1+) 



1 

X (4-|(2 + 3)(7+l)|6+) ' (125) 
((7+| (3 + 4) |2+) (54) + (7+| 6 |5 + ) (2 4)) 4 

Cl36 " 



(2 3) (3 4) (4 5) (5 6) [7 1] (1 + | (2 + 3) |4+) (7 + | (5 + 6) |4+> 



X (6-|(7+l)(2 + 3)|4+)(2-|(7+l)(5 + 6)|4 + ) ' (126) 
c = <2 5> 4 [6 7] 3 

147 ' S671 (23)(34)(45)[71](l + |(6 + 7)|5 + )(6+|(7+l)|2+) ' 1 1 

(3 + | (1 + 2) |5 + ) 4 

C236 = si 23 [12] [2 3] (4 5) (5 6) (6 7) (3+| (1 + 2) |7+) (1+| (2 + 3) |4+) ' (128) 

(5l(3 + 4)(6 + 7)|5+) 4 

C247 " 



[1 2] (3 4) (4 5) (5 6) (6 7) (2+| (3 + 4) |5+) (1 + | (6 + 7) |5+) 

X (5-| (6 + 7)(1 + 2) |3+) (5- 1 (3 + 4)(1 + 2) |7+) ' (129) 

(12) 3 [4 6] 4 

256 S4 56 (23)[45][56](71)(4+|(5 + 6)|7 + )(6 + |(4 + 5)|3 + ) ' 1 ) 
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(12) 3 (3 5) 4 [6 7f 

C257 



(2 3) (3 4) (4 5) (6 + | (4 + 5) |3+) (7+| (1 + 2) |3+) 



1 

X (3-| (4 + 5)(6 + 7) |1+) (3- 1 (1 + 2) (6 + 7) |5+) ' (131) 
(12) 3 (4+|(6 + 7)|1+) 4 

C357 



(2 3) [4 5] (6 7) (7 1) (4+| (2 + 3) |1+) (5+| (6 + 7) |1+) 

X (3~| (4 + 5)(6 + 7) |1+) (6- 1 (4 + 5) (2 + 3) |1+) ' (132) 
. = (12) 3 [6 7] 3 

367 ~ S567 (23) (34) [56] (5 + |(6 + 7)|l+) (7 + |(5 + 6)|4+) ' 1 ) 

Note that cd and C256 are flip-symmetric. 

The quantities (122)— (133) obey the linear relations (plus their flips), 

CA + C357 + C 367 = C B + C 25 6 + C 2 57, (134) 
C A + C C - C U 7 - C236 + C 3 57 + C 3 67 = [ca + C C ~ C U 7 ~ C236 + C 3 57 + C 36 7] | flip , (135) 

cc + C257 = [cb - c 23 e] I flip + c M 7 + C247, (136) 

cd + ci 36 + c 367 = [c c + c 357 ] I flip + C147, (137) 

ca - ci 36 - c 236 = [ca - ci 36 - c 236 ]|flip • (138) 
The tree amplitude can be written in a manifestly flip-symmetric form as 

CC = C 247 + C 2 56 + C B + C147 + [c B + C147] Iflip • (139) 

The remaining independent box coefficients are given by 



C123 


= [Ca] flip + Ci 3 6 + C 23 6 + C 2 56, 


(140) 


C127 


= C C + C 25 7 + [C236 + C367] Iflip , 


(141) 


Cl 3 4 


= C D + C147, 


(142) 


Cl57 


= c c + [c 367 ] | fli p, 


(143) 


C235 


= C A + C 2 56, 


(144) 


C237 


= c B + c 367 , 


(145) 


c 267 


= Cb + C2 3 6, 


(146) 


C345 


= C C + Ci 47 + [C147] |flip + C 357 , 


(147) 


C356 


= c A + c 236 , 


(148) 


c 457 


= Cc + c M7 , 


(149) 


C567 


= CA + C 236 + C 357 + C 367 . 


(150) 
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The final 16 box coefficients are obtained by the flip symmetry via eq. (71), where 

i = (4 - i) mod 7, (151) 
and again ijk should be written in ascending order. 

D. ^= 4 (l-,2+,3-,4+,5-,6+,7+) 

For this partial amplitude we define the "flip" operation, 

-X" | flip = —X(l ^5,2^4,6^7), (152) 
which is a symmetry of the tree amplitude, 

4 ree (l-, 2+ 3", 4 + , 5", 6+, 7+) = 4 rcc (l-, 2+ 3", 4+, 5" 6+ 7+) , (153) 

flip 

and of the one-loop amplitude. The box coefficients c^fc are given in terms of the following 
quantities: 

(13) 4 (7+|(4 + 6) |5+) 4 

° A ~ s 123 s i56 (1 2) (2 3) (4 5) (5 6) (7+ 1 (5 + 6) |4+> (7+ 1 (1 + 2) |3+> (6" | (4 + 5) (2 + 3) 1 1+) ' 

(154) 

r (13) 4 (4+| (6 + 7) |5+) 4 

B " 3123*567 (1 2) (2 3) (5 6) (6 7) (4+| (5 + 6) |7+> (4+| (2 + 3) |l+> (3"| (1 + 2)(6 + 7) |5+> ' 

(155) 

r = (l-|(6 + 7)(2 + 4)|3 + ) 4 

° " ^234*671 (2 3) (3 4) (6 7) (7 1) (5+| (3 + 4) |2+> (5+| (6 + 7) |1+) (4"| (2 + 3) (7 + 1) |6 + ) ' 

(156) 

-(5-|(4 + 6)(7 + 2)|l+) 4 

S456S712 (1 2) (4 5) (5 6) (7 1) (3 + | (4 + 5) |6+> <3+| (1 + 2) |7+) (4"| (5 + 6) (7 + 1) |2+> ' 

(157) 

((7+| (2 + 4) |3+) (54) + (7+| 6 |5+) (34)) 4 
(2 3) (3 4) (4 5) (5 6) [7 1] <1+| (2 + 3) |4+> (7 + | (5 + 6) |4+> 

1 

X (4- 1 (5 + 6)(7 + 1) |2+) (4- 1 (2 + 3)(7 + 1) |6+) ' (158) 
(2 + | (1 + 3) |5 + ) 4 

C236 = Sl23 [1 2] [2 3] (4 5) (5 6) (6 7) (3+| (1 + 2) |7+> (1+| (2 + 3) |4+> ' (159) 
((4+|(2 + 7)|l + )(5 7) + (4+|6|5 + )(17)) 4 

C246 



CD = 



Cl36 



(1 2) [3 4] (5 6) (6 7) (7 1) (4+| (5 + 6) |7 + ) <3+| (1 + 2) |7 + ) 



X (7-|(l + 2)(3 + 4)|5+)(7-|(5 + 6)(3 + 4)|2 + ) ' (160) 
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c = (13) 4 [4 6] 4 mi) 

256 S456 (1 2)(23)[45][56](71)(4 + |(5 + 6)|7 + )(6+|(4 + 5)|3 + ) ' 1 ) 



;i3) 4 (3 5) 4 [6 7] 3 



C257 = 



;i 2) (2 3) (3 4) (4 5) (6 + | (4 + 5) |3+) (7+| (1 + 2) |3+) 



1 

X (3-| (4 + 5)(6 + 7) |1+) (3- 1 (1 + 2) (6 + 7) |5+) ' (162) 

[2 4] 4 (15) 4 

0347 " s 23 4 [2 3] [3 4] (5 6) (6 7) (7 1) (2+| (3 + 4) |5 + ) (4+| (2 + 3) |1+) ' (163) 



C357 = 



;i3) 4 (4 + |(6 + 7)|l+) 4 



(1 2) (2 3) [4 5] (6 7) (7 1) (4+| (2 + 3) |1+) (5 + | (6 + 7) |1+) 



X (3-|(4 + 5)(6 + 7)|l+)(6-|(4 + 5)(2 + 3)|l + ) ' (164) 
r = (13) 4 [6 7] 3 

367 S567 (1 2)(23)(34)[56](5 + |(6 + 7)|l + )(7 + |(5 + 6)|4 + ) ' 1 ) 

Note that cp, C257 and c 347 are flip symmetric. 

The linear relations obeyed by the quantities (154)-(165) are 

CA + C357 + C 36 7 = C B + C 256 + C 25 7, (166) 

C C + C D + Ci 36 + C 367 — C B — C 2 36 + C 2 56 + [cb ~ C236 + C 2 5g] (flip + C 2 57 + C 3 47, (167) 

C C + Cl36 + C 36 7 = [CC + C 136 + C 367 ] | flip , (168) 

CB - C 2 36 - C 2 46 + C 3 47 = [cc - C357] | flip , (169) 

plus the equations related by the flip symmetry. The simplest, manifestly flip-symmetric 
representation we have been able to find for the tree amplitude is 

Aj C ° = C 257 + C347 + CB + C 256 + [CB + C256] | flip • (170) 

The remaining independent box coefficients are given by 



C123 = 


CD + C136 + C 2 36 + C 2 56, 


(171) 


Cl26 = 


CD + C256; 


(172) 


Cl34 = 


CC + [C367]|flip , 


(173) 


Cl37 = 


CC + C347, 


(174) 


Cl67 = 


CC + C136 + C347 + C 3 67, 


(175) 


C234 = 


C B + C347 + C 36 7 + [C357] flip , 


(176) 


C235 = 


CA + C256; 


(177) 
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C 2 37 = C B + C 367 , (178) 

C 2 67 = C B + C 236 , (179) 

C 3 56 = C A + C 23 6, (180) 

C567 = C A + C 23 6 + C 3 57 + C 3 67- (181) 

The final 16 box coefficients are obtained by the flip symmetry (71), where 

l = (7 - i) mod 7, (182) 
and again ijk should be written in ascending order. 



V. CONSISTENCY OF THE RESULTS 



In this section we describe various consistency checks that we performed on the ampli- 
tudes. Amplitudes, in general, must satisfy a stringent set of constraints on their analytic 
properties. In particular, all kinematic poles in the amplitudes must correspond to the phys- 
ical propagation of particles. Such constraints are so restrictive that they can be used, for 
example, to construct ansatze for infinite sequences of all-plus-helicity amplitudes in gauge 
and gravity theories [42, 52]. The poles may be divided into three categories: collinear, 
multi-particle and spurious. We consider each of these categories in turn. 



A. Collinear Behavior 

An important constraint on the amplitudes arises from the region in phase space where 
the momenta of two legs a and b become collinear. In the collinear region, k a — > zkp, 
k b — > (1 — z)kp, where kp is the momentum of the quasi-on-shell intermediate state P, with 
helicity A. In this limit, massless color-ordered tree amplitudes behave as 

AT ^ 5ZSplit tr ^ c (z,a Aa ,6 A6 )A t n ^ c 1 (...(a + 6) A ...), (183) 

A=± 

where Split ^ e are tree- level splitting amplitudes [18]. At one loop, the generalization is, 

A=± 

+ Split ^ oop (^, a Xa , b Xb ) ATi(- ■ ■ (a + b) x . . .)) , (184) 
31 



where the Split ^ oop are one- loop splitting amplitudes, which are tabulated in the second 
appendix of ref. [9]. This reference also contains a discussion of the behavior of the collinear 
limits of one-loop amplitudes and integral functions. 

We have verified numerically that all helicity amplitudes presented in section IV have the 
proper collinear behavior in all channels, as dictated by eq. (184). This, by itself, provides 
a rather powerful check on their form. 



B. Mult i- Particle Factorization 



A related constraint is that of multi-particle factorization. This factorization plays a 
central role in the tree-level CSW construction. Presumably, it will also play an important 
role in extending the CSW construction to loop level. These properties are not as well- 
discussed in the literature as the collinear properties, so we present the key points here, and 
include an example in appendix B. 

At tree level, multi-particle factorization is reasonably straightforward. For (fcj + ki + \ + 
• • • + /cj +r -i) 2 = K 2 — > (with r > 2), the amplitude behaves as 

^tree > Af^\(kj, . . . , kj +r -i, K x ) A^ r+1 ( (— K) A , fcj+ r , . . . , , (185) 

A=± 

where A is the helicity of the intermediate state with momentum K. 

At loop level, in infrared-divergent gauge theories, multi-particle factorization is more 
subtle. For such theories, loop amplitudes do not factorize in any naive sense, due to 
the emission of soft gluons which induce non-trivial logarithmic corrections to factorization 
formulae. These logarithmic corrections do however obey universal formula? [53] because of 
their close relation to the universal infrared divergences [50]. More explicitly, for K 2 — > 
the factorization properties for one-loop amplitudes are described by [53], 



,1-loop k 2 ->o 



A=± 



A r+1 .±(ki, . . . , ki +r -±, K ) — - A n _ r+l {{ K) , k i+r , . . . , 



K 2 



+ A^\(ki, . . . , fcj +r _i, K X ) ^n-r+l;l(( — ^0 A ) ^i+r, • • • , (186) 

+ A^Kki, . . . , fcj+ r _i, K x ) A^^ r+1 {{— K) A ; ki +r , . . . , cyT n {K 2 ] k±, . . . , k n ) 

where the one-loop factorization function T n is independent of helicities. This formula is 
similar to the one for an amplitude which factorizes naively, depicted in fig. 3, except that 
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i + r - 1 




FIG. 3: A schematic depiction of multi-particle factorization at one loop. 

J-'n contains kinematic invariants with momenta from both sides of the propagator carrying 
momentum K. For example ln(— Sj-i^) is one such logarithm: is a momentum belonging 
to one of the factorized amplitudes on the right-hand side of eq. (186) and ki to the other 
amplitude. 

In general, the factorization function is composed of 'factorizing' and 'non-factorizing' 
components, 

-77 Tfact i Tnon-fact /i 
•>-n = J- n +S n ■ (ISO 

The factorizing contributions are easily obtained by computing bubble Feynman diagrams. 
For the case of Af = 4 supersymmetric gauge theory the factorizing contributions vanish, 
leaving only the 'non-factorizing' contributions. 

These 'non-factorizing' contributions are linked to the infrared divergences, as shown 
in ref. [53]. A constructive proof of the universality of these non-factorizing contributions 
was also given in that reference, as well as explicit formulae for determining their values in 
any theory. All non-factorizing contributions are determined from the infrared divergences 
present in a given process. The infrared divergences in the J\f = 4 theory at one loop are 

n / 2 

_cr / H 

e pole 



<i =4 (l,2,...,n) 



__£ 



A 



tree , 



1,2, 



,n) . 



(188) 



From Table 2 of ref. [53], the appearance of the singularities, 



cr 

c 2 



cr 

c 2 



— l,i+r 



A 



tree 



(189) 



0, the factorization function is 
2 



implies that in the factorization channel K = Sj...(j+ r -i) 

A {K 2 \ hi, . . . , k n ) = 2(fi 2 ) e (^F^\. i+l + -F^. r ^ r _ 1;i+f , +1 + — (— Sj.. .( i+r _i)) e j 
— 2(^ 2 ) e ^F 2m?i (sj_i i i, Sj...( i+r _i), S( i+ i)...( i+r ._i), S( i+r )...(j_2)) 
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_|_ _p 2mfe ( Si+r _ 1 ^ +r? s i ___^ +J1 _ 1 j ) S(j +r+ i)...(j_i), Sj...(j+ r -2)) 




(190) 



where we use the notation from ref. [53] on the first line. On the subsequent lines we use 
the notation of appendix A for the integral functions. 

We have checked that all amplitudes in section IV satisfy the correct multi-particle fac- 
torizations dictated by eq. (186), with the factorization function (190). In appendix B, we 
present an example of a multi-particle factorization limit. All the remaining multi-particle 
factorization limits, for any of the seven-point helicity amplitudes given in section IV, are 
similar. 

C. Cancellation of spurious singularities 

In addition to the collinear and multi-particle poles discussed in the previous two subsec- 
tions, the coefficients also contain 'planar' and 'cubic' singularities as defined in section III. 
Unlike the collinear and multi-particle singularities, these cannot be singularities of the 
whole amplitude. Factorization (even in massless gauge theories) does not allow them. We 
therefore expect them to cancel between different integral functions. Their origin in the 
integral reductions also suggests this. 

This is indeed what one finds. For example, in A^= 4 (l~, 2", 3", 4+ 5 + , 6+ 7+), there 
are eight coefficients that contain a denominator factor of (1 + | (2 + 3) |4 + ). As noted in 
section III, these coefficients — ci 36 , ci 67 , C236, C267, C346, C356, c^-j, and — will therefore 
be separately singular when fc 2 + k 3 = a\k\ + a 4 /c 4 . One can nonetheless verify that these 
singularities cancel in the amplitude as a whole. Indeed, the cancellation occurs separately 
within two linear combinations of terms in the amplitude, 



While a subset of terms in the whole amplitude will allow the cancellation of any given 
singularity — typically four terms for the planar singularities, and five for the cubic ones - 
there does not appear to be any subset of terms in which all spurious singularities cancel. 

We have checked numerically that all spurious singularities do indeed cancel in the entire 
amplitude, for each of the four NMHV helicity configurations. This cancellation is a strong 



c 136 B(l, 3, 6) + c 236 B(2, 3, 6) + c U6 B(3, 4, 6) + c 356 B(3, 5, 6) , 



and 



c 167 B(l, 6, 7) + c 267 B(2, 6, 7) + c 467j B(4, 6, 7) + c 567 B(5, 6, 7) . 



(191) 
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check on the amplitude, precisely because it involves the analytic behavior of the integral 
functions in a non-trivial way. 



VI. AN ALL-MULTIPLICITY COEFFICIENT 

The computation of amplitudes beyond a fixed number of external legs can provide ad- 
ditional information and clues to general structure. In particular, as we shall discuss in sec- 
tion VII, a discussion of the twistor-space structure of the seven-point amplitude requires an 
accounting for the 'holomorphic anomaly' pointed out by Cachazo, Svrcek, and Witten [37]. 
Some aspects of the structure of an all-n amplitude also require such an accounting; but we 
will be able to describe most aspects without reference to the 'anomaly'. 

We limit ourselves here to confirming the coefficient presented in ref. [40]; presenting one 
new set of coefficients, of a class of three-mass box functions; and discussing a number of 
vanishing coefficients. 

Consider the adjacent-minus amplitude A n -i(l~, 2~ , 3~, 4 + , . . . , n + ). The simplest of the 
cuts is the cut in the S123 channel, represented by the product of tree amplitudes, 

Cm = i I", 2" 3^, tt) Af* ((-£4)-, 4 + , 5 + , . . . n + , £~) , (192) 



integrated over phase space. This cut is particularly simple because only MHV (or MHV) 
amplitudes participate, and only gluons can cross the cut (so the result is the same in J\f = 4 
or pure Af = 1 super- Yang-Mills theory, or QCD). Inserting eq. (7) twice, the cut becomes 

r = *( s i23) 3 1 nq ^ 

123 [12][23]<45)<56>...((n-l)n> [l x 1] [3 £4} (£4 4) (nl,) ' 1 ) 

We denote the running loop momenta by £f, that is, £ n = £\ + k n , £±, £2 = £\ — ki, £3 = £1 — 
ki — k 2 , £4 = £i — k 1 — k 2 — k 3 . To obtain eq. (193) we have used (£\ £4) [£ 1 £4] = [£\—£4) 2 = Si23- 
We can also clear out the loop-momentum-dependent spinor-products from the denominator, 
in favor of standard propagators, 

Cl23 = ^1! (3-| ft|4-)(l-| hW) 

123 [12][2 3](4 5>(5 6>...((n-l)n> Qfiftft ' 1 J 

It turns out that the algebraic steps for integrating the cut are essentially identical to those 
needed in the six-point case [10]. One can actually do a bit better, and avoid any integration, 
by multiplying and dividing by the spinor strings, 

(n~\ (1 + 2) |3-)(4-| (2 + 3) |l") , (195) 
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whose appearance is motivated by the form (22) of the denominators encountered in reducing 
non- adjacent two- mass pentagon integrals. One obtains, 

r = ^(g!23) 3 

123 [1 2] [2 3] (4 5) (5 6) • • • <(ra - 1) n) (1+| (2 + 3) |4+> (3+| (1 + 2) |n + ) 
i(l + 75 )(l + 2)3f 4 4(2 + 3)lf 1 nl 



tr 
x — 



p2 p2 p2 p2 



(196) 



The trace of Dirac 7 matrices, with a chiral projection inserted, is easily evaluated. The 75 
terms are proportional to the Levi-Civita tensor and integrate to zero. One can algebraically 
reduce the non-Levi-Civita terms to a combination of four cut box integrals, 



A 1 - loop (l-,2-,3-,4 + ,...,n 
= icr 



123 cut 

,3 



S 123 



[1 2] [2 3] (4 5) (5 6) • • • ((n - 1) n) (1 + | (2 + 3) |4+) (3+| (1 + 2) \n + ) 

s nl, s 123) s 23; s 4...(n-l)) 

+ F 2m/l (s 3 4,s 123 ,s 5 ... n ,s 12 ) 

+ F 2me { Sn ... 3 , S1...4, Si23, S 5 ...(n-1))) , (197) 

where the integral functions F are given in appendix A. This result reduces to eqs. (45) 
and (51) for n = 7. It also confirms Cachazo's [40] very elegant evaluation of this coefficient, 
obtained by exploiting the holomorphic anomaly [37]. 

Now consider the coefficient of a class of three-mass box functions in the same amplitude. 
In this case NMHV tree amplitudes appear in the cuts. We used the simplified form of the 
NMH V amplitudes appropriate for use in the cuts, obtained in ref . [30] . The cut starts out 
with the structure of a tensor octagon integral (independent of n), and can be cleanly reduced 
to a sum of scalar box integrals. Here we present the coefficient of the set of three-mass box 
integrals, 

-F 3m (.S2...(ci-l) ? S( C2 +1)...2, S 3 ...( ci _i), S C1 ... C2 , S( C2 +l)...l) , (198) 

shown in fig. 4. This coefficient can be written as a single term, 

((12) (2 3) Sc ,„ C2 ) 4 ((ci-l)c!) (eg (c 2 + l)) 

r <12>(23>..>1> * Sc ,.. C2 ((ci-l)1^ c ,.. C2 ^ (c2+1) ... 1 |2+)( C r|^ c ,.. C 2^(c2 + i)...i|2 + > 

1 



x 



(C 2 -|^ C ,.. C2 ^ 3 ...( C1 -1) |2 + ) ((C 2 + 1)-|^ C ,.. C2 ^3...( C1 -1) |2 + ) ' 

(199) 
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FIG. 4: The class of three-mass box functions whose coefficient is given in eq. (199). 

We have checked that this set of coefficients correctly reproduces the appropriate box coef- 
ficients for n = 6 and 7. 

We defer a more detailed discussion of this computation, as well as results for other 
box coefficients, to a future publication. While the all-n NMHV loop amplitude contains 
everything from one- mass to three-mass boxes (four-mass boxes are absent), a generic term 
contains either a three-mass box, or an easy two-mass box. The latter coefficients are more 
complicated in structure, while the former are simpler. 

Using the triple cut shown in fig. 1, it is easy to see that the coefficient of any box 
integral which contains two adjacent clusters of legs, each having only positive-helicity gluons, 
must vanish. This statement is true for any amplitude in any cut-construct ible theory. 
The reasoning is as follows: Consider a triple cut where each of the two positive-helicity 
clusters is identified with the set of external legs emerging from a blob in fig. 1. Then the 
vanishing of the tree amplitudes A^ cc (±+H — • +) implies that the particle crossing the cut 
line between these two blobs must have (— ) helicity on both sides of the cut, which is not 
a valid helicity assignment. A corollary is that any box integral with one all-plus cluster 
adjacent to two consecutive massless plus-helicity legs must have a vanishing coefficient. 

A simple consequence of this result for the seven-point amplitude ( h+++) is that 

C235 = C237 = C257 = C357 = 0, as given in eq. (50). Clearly a large number of box coefficients 
for the all-n amplitude ( H h+) vanish by these rules. 
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VII. TWISTOR-SPACE PROPERTIES 



The target space for Witten's candidate topological string theory is CP 3 ' 4 , otherwise 
called projective (super-)twistor space. Points in twistor space correspond to null momenta 
or equivalently to light cones in space-time. The correspondence is specified by a 'half- 
Fourier' transform. More precisely, if we represent a null momentum by the tensor product 
of a spinor A" and a conjugate spinor A a , then twistor quantities are obtained by Fourier- 
transforming with respect to all the A a . 

Amplitudes in twistor space, as it turns out, have rather simple properties. At tree level, 
they are non-vanishing only on certain curves. This implies that they contain factors of delta 
functions (or derivatives thereof) whose arguments are the characteristic equations for the 
curves. The coefficients of the delta functions, however, have been quite difficult to calculate 
directly. 

As Witten pointed out in his original paper [1], however, we do not need the twistor- 
space amplitudes in order to establish the structure of the delta functions they contain. In 
momentum space, the Fourier transform turns the polynomials into differential operators 
(polynomial in the Aj, and derivatives with respect to the Aj), which will annihilate the 
amplitude. One particularly useful building block for these differential operators is the line 
annihilation operator, expressing the condition that three points in twistor space lie on a 
common 'line' or CP 1 . If the coordinates of the three points, labeled i, j, k, are Z\ = (A", /x"), 
etc., then the appropriate condition is 

e IJKL z!Z J j Z* = Q, (200) 

for all choices of L. Choosing L — a, and translating this equation back to momentum space 
using the identification fi a <-> —id/dXa, we obtain the operator, 

Fijk = (ij) J- + (jk) ^- + (hi) J- . (201) 
o\k oXi oXj 

Two important sufficient conditions for to annihilate an expression, i. e. for it to have 
support only when k lie on a line in twistor space, are [1] 

1. The expression is completely independent of Aj, Xj, and A&, or 

2. Aj, Xj, Afc appear only via a sum of momenta containing them, of the form 

P ad = (••• + 1* + k, + h + • • -) ad = ■■■ + XfXf + A, a Aj + XfXt + ■■■ . (202) 
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The first condition is obvious from the definition (201); the second holds because of the 
Schouten identity, 

(ij)\k+{jk)\ i + (ki)\ j = 0. (203) 

The tree-level MHV amplitude, for example, is annihilated by because it is indepen- 
dent of the Aj. Any possible delta functions vanish for generic momenta, because they take 
the form 5({ij)). At one loop, Cachazo, Svrcek, and Witten [37] pointed out that such delta 
functions, arising from the spinor analog of the fact that dz(l/z) ^ 0, do arise. They must 
be taken into account for a proper analysis of the twistor-space structure of amplitudes. 

We will not compute the relevant 'holomorphic anomaly' terms for the amplitudes in 
this paper, and so we will not be able to fully exhibit their twistor-space structure. While 
the 'anomaly' terms enter into the action of the differential operators on the box integrals, 
their action on the coefficients is unaffected by it. The properties of the coefficient are also 
important, and we discuss them. 

In addition to the line operator we will employ the planar operator [1], 

K m = e IJKL Z{Z J .ZK = (ij) e^-^-^ ± [5 permutations] , (204) 

o\ a k dxl 

whose vanishing implies that four points lie in a plane in twistor space. 

In the case of the seven-point amplitude, all line operators F and planar operators K 
are affected by the 'anomaly'. For higher-point amplitudes, however, there are classes of 
unaffected operators. We will study their action on the term in the amplitude discussed in 
section VI. We shall show it is consistent with the simple twistor-space structure expected 
from a generalization of the BST calculation. 

A. Properties of seven-point box coefficients 

Before describing the twistor-space structure of the coefficients from section IV, we 
list some of their other (related) properties here. 

• Each box coefficient can be written as the sum of a small number of 'terms', or 'one- 
term coefficients': namely, the auxiliary quantities ca,cb,cc,cd,ce (as needed), plus 
the Cijk that are given directly in section IV in terms of spinor strings. 

• The sum contains at most four terms. In the case of A^= 4 (l~, 2~, 3~, 4 + , 5 + , 6 + , 7 + ), 
it has at most three terms. 
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Only the coefficients of the one-mass box integrals, (with indices mod 7), 

require the use of three or four terms. 

The coefficients of the hard two-mass box integrals, 0^+1^+3 or Q_2,j,j+i, require at 
most two terms. 

The coefficients of the easy two-mass box integrals, 0^+1^+4, and of the three-mass 
box integrals, 0^+2,1+4, are all one-term expressions. 

Each term resembles a product of off-shell MHV vertices [23] , in terms of the types of 
denominators that occur: multi-particle invariants but not full two-particle 

poles s^j+i; spinor products continued 'off-shell', like (aB*) = (a~\ (b + c) \d~), where 
k-B = k b + k c , and d temporarily plays the role of the arbitrary reference spinor in the 
CSW formalism (see also ref. [30]). The choice of d varies from term to term, however. 
The longer strings of the form (a~\ {b + c)(d + e) |/ + ) could correspond to continuing 
both spinors off-shell. 

The resemblance extends to the numerator factors, which always appear raised to the 
third or fourth power, like the factor of (ij) in the numerator of the Parke- Taylor 
amplitudes [7] . Probably they should always be considered raised to the fourth power, 
because whenever such a factor appears raised to the third power, it always seems to be 
of the type that 'could have appeared' in the denominator, such as (i, i + 1), [i, i + l], 
Si,i+i,i+2, or even the more complicated spinor strings which are 'square roots' of the 
integral reduction factors discussed in section III. Other numerator factors which 
cannot occur in denominators, such as (i, i + 2), or (a~ \ (b + c) \d~) where a, b, c, d are 
not cyclicly consecutive, always appear raised to the fourth power, if they appear in 
the numerator at all. 

Perhaps the most interesting of the one-term coefficients are those exemplified by C135 
in Aj[^ 4 (l~, 2~, 3 + , 4~, 5 + , 6 + , 7 + ). Here a new structure appears in the numerator, 
unlike any of the denominator factors. The unexpected form made the simplified 
form of this coefficient more difficult to guess. In the end, we deduced it from its 
collinear behavior in particular channels. Using the Schouten identity and momentum 
conservation, one can show that it has the following collinear limits, 

(3 + |(6 + 7)|l + )(46) + (3 + |5|4+)(16) -> + <3+| (6 + 7) |l + ) (46) for 4||5, (205) 
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-> - (3+ 1 (2 + 4) 1 1+) (4 6) for 5||6, (206) 
-> - <3 + | (1 + 2) |4+) (16) for 6||7, (207) 
-> + (3+ 1 (5 + 6) |4 + ) (1 6) for 7| 1 1. (208) 

The most general twistor-space property of the box coefficients is that all points lie in a 
plane. That is, K mnpq for every choice of m, n, p, q annihilates every one-term coefficient, and 
hence, by linearity, it annihilates every box coefficient c^. This property is not so easy to 
see analytically, but it is straightforward to verify numerically by evaluating the expressions 
K mn pq c ijk at random points in 7-particle phase space. 

On the other hand, the way the 7 coplanar points get distributed into lines varies from 
term to term, yet it does so in a very systematic way. Also, the correct behavior can be 
determined by inspection. From the amplitude A^f 4 (l", 2", 3 + , 4+ 5 + , 6 + , 7 + ), consider, 
for example: 

• the easy two- mass box coefficient C125 given in eq. (78). Here legs 3,4,5,6 only appear 
as Xi spinors, not Aj spinors. So legs 3,4,5,6 must lie on a line for C125. On the other 
hand, legs 7,1,2 appear sometimes as A« spinors, as in [1 2] and [71], so they should be 
placed at generic points in the plane, as shown in fig. 5a. 

• the hard two-mass box coefficient C237 given in eq. (85). Here legs 1,2,3 either appear 
as Aj spinors, or through the momentum sum k\ + k 2 + k 3 . The latter appearance 
is easy to see in s 12 3 = (&i + k 2 + k 3 ) 2 , but it is also true for (4 + | (2 + 3) |1 + ) = 
(4+| 7 m|i+) (k 1 + k 2 + k 3 Y, since ^|1+) = 0; and similarly for (3~| (l + 2)(6 + 7) |5+). 
Thus legs 1,2,3 lie on a line. Likewise, legs 5,6,7 lie on a different line. Leg 4 appears 
as A4, so it sits in the plane, away from the 2 lines, as shown in fig. 5b. 

• the three-mass box coefficient C135 given in eq. (79). From the denominator factors 
alone, one would conclude, as in the previous case that legs 4,5,6 lie in a line, and legs 
6,7,1 lie in a line. Thus the two lines intersect at a common point, leg 6. The form 
of the numerator in eq. (79) makes manifest the proper 6,7,1 behavior. However, it is 
not apparent that it allows 4,5,6 to lie in a line until one rewrites 

<3+| (6 + 7) |l + ) (46) + <3 + | 5 |4+> (16) = <3 + | (5 + 6) |4+) (16) + <3+| 7|l+) (46) . 

(209) 
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(a) (b) (c) 

FIG. 5: Examples of twistor-space configurations for single-term box coefficients in the helicity 
amplitude A^f 4 (1~, 2", 3+, 4", 5+, 6+, 7+). In every case, all the points lie in a plane, (a) the 
easy two-mass box coefficient C125, (b) the hard two-mass box coefficient C237, (c) the three-mass 
box coefficient C135. 

Numerical investigation reveals that the pictures shown in fig. 5 are accurate, in that there 
are no 'hidden' twistor-space co-linearities. Note also that each configuration in fig. 5 can 
be described as follows: all seven points lie in a plane, distributed among three lines, and 
two of the lines can be chosen to intersect at one of the seven points. 

In fact, the configurations shown in fig. 5, after permuting the external leg labels around, 
are the complete set of twistor-space configurations for all the one-term coefficients for all 
four helicity amplitudes. Furthermore, 

• The three-mass box coefficients (always one term) always have the configuration 
of fig. 5c. The massless leg in the corresponding integral always lies at the intersection 
of the two lines. 

• The easy two-mass box coefficients (always one term) always have the configuration 
of fig. 5a. The three points lying off the line in fig. 5a are those belonging to the 
three-leg cluster forming one of the external masses in the integral. 

• The hard two-mass box coefficients are the sums of (at most two) terms with the 
configurations of fig. 5a and fig. 5b. In the fig. 5a term, the three points lying off the 
line are those belonging to the two- leg cluster and the massless leg adjacent to it. In 
the fig. 5b term, the point lying off of both lines is the massless leg adjacent to the 
three-leg cluster. 

• The one-mass box coefficients generally involve all three types of configurations in 
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FIG. 6: The twistor-space configuration for the class of three-mass-box coefficients in the helicity 
amplitude ( h+ • • • ++) discussed in the text. All points lie in a plane. 

fig. 5. However, for all helicity amplitudes the full coefficients do obey ^712-^456^23 = 
and i 7 'i23-^456 c i23 = 0, plus all cyclic permutations thereof. 

The coefficients of the box integral functions in the six-point NMHV amplitude [10] have 
an analogous structure, only simpler because there are fewer points: All six points lie in 
a plane in twistor space, and for each single term in a box coefficient, three of the points 
always lie on a line. 

B. Properties of an all-n box coefficient 

The result for the class of three-mass-box coefficients for the all-n NMHV amplitude 

( h+ • • • ++) is given in eq. (199). We now analyze its twistor-space structure. First 

observe that all points in the set {ci,ci + 1, ...,02} (where the indices are understood 
cyclicly modn) lie on a line, because they only appear through the momentum sum K Cl ___ C2 
or through A C1 or A C2 . The points {c 2 + 1, c 2 + 2, . . . , 1, 2} also lie on a line, for a similar 
reason. (Note that W( C2 +i)...i\^ + ) = -$'( C2 +i)...2|2 + ), and that leg 2 otherwise only appears as 
A2.) Finally, the points {2, 3, . . . , c\ — 1} lie on a line. Numerical investigation for a number 
of randomly- chosen box coefficients for n = 8, 9, 10 confirms that all n points lie in a plane. 
This twistor-space structure is summarized in fig. 6. 

This structure generalizes the behavior of the seven-point 'one-term' coefficients, which 
was illustrated in fig. 5. Again a 'one-term' expression can be described in twistor space as 
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having all points lying on three lines within a single plane; and one of the intersections of 
the lines always contains one of the n points. 

In the full amplitude, the coefficient (199) appears multiplied by the three-mass-box 
integral function -F 3m (s 2 ...( Cl -i), S( C2 +i)...2, s 3 ...( Cl -i), s ci...c 2 , s C2 +i...i) in e Q- ( 19 8). The action 
of the F and K twistor-space operators on F 3m is affected by the holomorphic anomaly 
The behavior of leg 2 in particular is affected, because of singularities in the integral when 
the loop momentum becomes collinear with k 2 . At the level of the F operator, none of 
the other legs are affected by this issue, since they are safely buried into massive clusters. 
So, the lines in fig. 6, excluding the point 2 and ignoring their intersections, also describe 
the twistor-space collinear behavior of the box function multiplying the coefficient. (This 
property is also manifest from the form of the arguments of F 3m .) On the other hand, the 
co-planarity of (the finite parts of) this integral is not apparent, and perhaps doubtful, even 
taking into account the holomorphic anomaly. 

VIII. CONCLUSIONS 

In this paper we computed all the next-to-MHV one-loop seven-gluon amplitudes in 
M = 4 super- Yang-Mills theory, as well as a class of contributions to an n-gluon next-to- 
MHV amplitude. These results provide a useful guide to the spinor-product structure and 
(related) twistor-space properties of higher-point M = 4 amplitudes. The coefficients of the 
box integral functions appearing in the amplitudes can be written as simple terms built out 
of spinor strings, or sums of just a few such terms. Each term exhibits a simple twistor-space 
structure: all points lie on three lines confined to a single plane, and two of the lines always 
intersect at one of the n points. 

The co-planarity of the NMHV box coefficients is a very intriguing result. This complete 
co-planarity may be demonstrated for general one-loop NMHV amplitudes along the same 
lines used by Cachazo to demonstrate a certain degree of co-linearity [40]. Using the co- 
linearity of MHV tree amplitudes, Cachazo showed that F^ annihilates a box coefficient if 
the points i, j, k belong to the MHV side of a 'nontrivial' cut, one for which the corresponding 
box function has an imaginary part. If the coefficient were not annihilated, logarithms would 
remain in the cut after applying Fijk, in contradiction to the result obtain by integrating 
the holomorphic anomaly in the cut [39, 40]. (Because the anomaly delta function freezes 
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the phase-space integration, the action of can give only a rational function or zero.) 
The same logic can be used to determine all and F^ that annihilate a given box 

coefficient. On one side of a cut the tree amplitude is MHV and therefore supported on 
a line. On the other side of the cut the tree amplitude is NMHV so it is supported on 
pairs of intersecting lines. From the lack of logarithms in the integrated anomaly, this 
implies that the coefficient of any integral with a cut in a given channel satisfies the same 
coplanarity and collinearity requirements satisfied by the tree amplitudes on either side of 
the cut. The complete co-planarity of a box coefficient can be demonstrated by combining 
the information from overlapping cuts of the associated integral. This simple argument 
works for all coefficients, except for those for one-mass boxes, or for easy two-mass boxes 
when one of the masses is composed of only two external legs. In these cases, however, one 
can appeal to the linear equations obtained from the infrared singularities, eq. (28), for two- 
and three-particle invariants, in order to argue (for n > 6) that these coefficients must also 
be co-planar. 

This twistor-space co-planarity is stronger than would seem to be implied by generic 
MHV loop diagrams. For NMHV one-loop amplitudes these diagrams have three vertices. 
The one-particle-irreducible MHV-diagrams (triangle diagrams) heuristically appear to be 
planar. (An MHV vertex is dual to a line in twistor space; a propagator in the loop is dual to 
an intersection between two lines [1, 23, 38].) On the other hand, the one-particle-reducible 
MHV-diagrams (bubbles with an additional MHV vertex attached to an external leg) do not 
appear to have to be planar. Apparently the sum is better behaved than individual MHV 
diagrams. 

The distribution of the seven points into lines in the plane, for the one-term coefficients in 
Ayf x =4 , as depicted in fig. 5, is also intriguing. This division is independent of the particular 
NMHV helicity configuration, and only depends on the type of integral being considered. 
This helicity-independence seems obscure from the point of view of MHV diagrams. 

The simplicity of the structure we have uncovered suggests that the computation of 
further multi-leg one-loop M = 4 amplitudes will be quite tractable, and provides strong 
motivation for additional investigations in this direction. 

The seven-point results presented here may also have some practical relevance. They 
can be used to benchmark purely numerical approaches to complicated multi-parton loop 
amplitudes in QCD, which are likely to become available in the near future. Such amplitudes 
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will contribute to a better understanding of cross sections for multi-jet production at hadron 
colliders. Multi-jet rates are important to understand for their own sake, but also because 
such processes form both irreducible and reducible backgrounds (the latter when jets fake 
leptons or photons, for example) to processes involving the electroweak interactions. We 
are optimistic that combining insights from unitarity and twistor space will provide new 
practical tools for computing experimentally-relevant loop amplitudes in massless gauge 
theories. 
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APPENDIX A: BOX INTEGRALS 



In this appendix we collect the dimensionally-regulated integral functions appearing in 
the M = 4 amplitudes; these integral functions were obtained from refs. [47, 55]. The reader 
is referred to these papers for further details of their derivation. Through O(e ), we have 
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where the ki denote massless momenta and the Ki massive momenta. The external momen- 
tum arguments K ly . . . ,K± are sums of external momenta ki from the n-point amplitude. 
The kinematic variables appearing in the integrals are 

s = (h + k 2 ) 2 , t = (k 2 + k 3 ) 2 , (A7) 

or with k relabeled as K for off-shell (massive) legs. The functions appearing in F 4 4m are 

P = \]l - 2Ai - 2A 2 + A 2 - 2AiA 2 + A 2 , (A8) 

and 
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We have rearranged the expressions for F 3m and F 2mh to make the poles in e more trans- 
parent. We have also corrected some signs in F 4m in ref. [9] and in the published version of 
ref. [47]. 

APPENDIX B: MULTI- PARTICLE FACTORIZATION EXAMPLE 



In this appendix we explicitly evaluate the multi-particle factorization of an NMHV 
seven-point amplitude in J\f = 4 super- Yang-Mills theory. As discussed in the text, one-loop 
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multi-particle factorization is described by a universal formula (186), even though naive 
factorization does not hold for infrared-divergent gauge theories. 

As an example, consider the amplitude A^ 4 (l~, 2 _ , 3~, 4 + , 5 + , 6 + , 7 + ). According to the 
general factorization formula, in the limit K 2 = s 2 34 — * 0, where K = k 2 + k 3 + A; 4 , we have 

A M=4 ^0 ^(2- 3- 4+> { _ K)+) -J- A *r(K-, 5\6\7\ 1-) 

S234 

+ Af e (2-, 3", 4 + , (-K) + ) — A«=\K~, 5+, 6+ 7+, 1") (Bl) 

S234 

+ AT C {2~, 3", 4+ (-K) + ) — Af c (ir-, 5 + , 6 + , 7+ 1") 

-§234 

X C r J^ =4 (s 2 34; fcl, • • • , ^7) • 

The subtlety in loop-level multi-particle factorization can be exposed by inspecting the 
infrared divergences on both sides of eq. (Bl). On the left-hand side, the divergences are [50], 

<f 4 (l,2,3,4,5,6,7) = -^W_^ ^(1,2,3,4,5,6,7). (B2) 

e pole e ^ \ - S j,j+l/ 

On the other hand, in the products of amplitudes appearing in the first two terms on the 

right-hand side of eq. (Bl), we have the singular terms 
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The mismatch between eq. (B2) and eq. (B3) must be absorbed into the factorization func- 
tion JF 7 appearing in the third term in eq. (Bl). This mismatch implies, after expanding in 
e, that the factorization function must contain the logarithmic terms 

ln (- fi i2) + ln(-s 45 ) ^ 

Neither of these terms allows for a separation of momenta {k 2 , k 3 , ^4} from {k 5 , fc 6 , kj, ki}, 
which one might have expected from a naive interpretation of eq. (Bl), or from fig. 3 with 
% — 2 and r = 3. 

The mismatched infrared singularities in eq. (B4) encode the factorization functions. 
From eq. (190), the infrared singularities (B4) imply that the factorization function in this 
channel is (recall F* =A = ^ on " fact ), 

^-^( S2U - kl ,...,k 7 ) 
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Now compare this prediction to the results for A^' 1 =4 (l~, 2~, 3~, 4 + , 5 + , 6 + , 7 + ) given in 
section IV. From that section, eight of the box functions have an s 2 34 kinematic pole in their 
coefficients. Thus, 



<f 4 (l-,2-,3-,4+,5+,6 + ,7^ 



S234^0 



— > c 134 5(1, 3, 4) + c 137 5(1, 3, 7) + c 167 B(l, 6, 7) + c 234 5(2, 3, 4) 

+ C345 5(3, 4, 5) + c 3 46 5(3, 4, 6) + C347 5(3, 4, 7) + c 46 7 5(4, 6, 7) . (B6) 

In this limit the coefficients appearing with an s 234 pole behave simply. Five of the coefficients 
give, 



{C134, c 234 , C345, c 346 , C347} -> Af°(2 ,3 ,4+, (-K)- 



while the remaining three give, 
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In the S234 factorization limit, an inspection of the explicit form of the integrals in 
appendix A reveals, after expanding in e, that the ones with an S234 pole in their coefficients 
behave as (see also table 5 of ref. [53]), 



B(2, 3, 4) = 5 lm (s 56 , s 67 , s 5 67) — > 5 lm (s 56 , s 67 , s 567 ) , 
5(3, 4, 5) = 5 lm (s 67 , s n , s 671 ) — > 5 lm (s 67 , s 71 , s en ) , 
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5(3, 4, 6) = F 2mh (s 71 , s 567 , S234 , s 56 ) — 5 lm (s 71 , s 567 , s 56 ) + ^(-s 2 3 4 ) 
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5(3, 4, 7) = F 2me (s 567 , s 671 , s 234 , s 67 ) — F lm (s 567 , s 671 , s 67 ) + ^(-s 2U )" , (B13) 
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The two remaining integrals, 



F m (si 2 , S 23 4, S34, S567), F m (S45, S 23 4, S671, S 23 ), 



(B15) 



do not have particularly simple properties as s 23 4 — > 0, although the polylogarithms do 
simplify to logarithms. 

Using eqs. (B7)-(B14), we may rewrite eq. (B6) as 



^f 4 (l-,2-,3-,4+,5+,6+,7 + ) 

— > c r (^YAT e (2-, 3", 4+, (-#)+) — A*~(K-, 5+ 6+, 7+, 1") 

S 2 34 



X 



F m (s 56 , S 67 , S 567 ) + F m (s 67 , S71, S 6 7l) + ^ m (s 6 71, S 5 6, «7l) 



(B16) 



+ F im (sn, S567, S 56 ) + F im ( S56 7, S 67 l, S 67 ) + 2F Um (s 23 , S 34 ) + ft 234 



where the remainder is 

5 



234 



■^y- S 234) ~ 2e2 



1 (-S234) e (-s 71 y 



-Li 2 (l- 



£7l\ 



1 (-S234) '(--§56)" 

2e 2 (-s 7 i)~ e " v s 7] 

,^2m /i / „ „ „ „ \ io 172m /i 



£ 2 - ^2 1 

+ 2F 2mh ( Sl2 , S234, S34, S567) + 2F 2mh (s 45 , S 234 , S 67 l, S 23 ) 



(B17) 



Comparing eq. (B16) to eq. (Bl), we may identify 



<f 4 (2~, 3", 4 + , (-K) + ) = 2 c r (/i 2 r 4 ree ^° m (^23, ^34) , 



(B18) 



in agreement with eq. (14), and 



<f 4 (^,5+,6+,7+,l-) = crV)^' 



2 \ e 4 tree 



F lm (s 56 , s 67 , s 567 ) + F lm (s 67 , s 7 i, s 67 i) 
+ F lm (s n , s 567 , s 56 ) + F lm (s 567 , s 67 i, s 67 ) 



+ F lm (s 6n , S 56 , S71) 



(B19) 
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in agreement with with eq. (15). The remainder may also be identified with the factorization 
function, 

•F 7 (s 234 ; h,..., hi) = (fi 2 Y K 2U . (B20) 

The last equation requires use of a dilogarithm identity and expansion through O(e ). 
Combining this result with eqs. (B18) and (B19), we thus find that the amplitude 
Ay^ =4 (l~, 2~, 3", 4 + , 5 + , 6 + , 7 + ) has precisely the expected factorization properties in the 
limit S234 — > 0. 

The other channels work similarly. For all the amplitudes in section IV, for each multi- 
particle limit the coefficients of precisely eight integrals have poles in the channel of interest. 
Five of the integrals belong with A^= 4 , one belongs with A^-f 4 , and the remaining two 
combine with the discontinuities from the other integrals to form the factorization function. 
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